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FIR FILTER AND SETTING METHOD OF COEFFICIENT THEREOF 

BACKGROUND OF THE INVENTION 

1. Fic^ld of the Invention 

The present invention relates to an FIR filter 
and setting method of coefficients of the FIR filter that 
are necessary for digital signal processing. 

2 . Description of the Related Art 

In the digital signal processing for picture 
and/or voice, filter processing is often used. Linear- 
phase FIR (Finite Impulse Response) filter is often 
utilized as the filter for digital signal processing of 
picture and/or voice in that the linear-phase FIR filter 
has the characteristics that its number of taps is finite 
and has linear-phase. 

Fig.l is a view illustrating circuit 
configuration of transversal type filter of the linear- 
phase FIR filter. 

The linear-phase FIR filter 1 , as illustrated in 
Fig.l, has (n - 1) delay units 2-1 to 2-(n-l) that 
constitute a shift register connected to an input 
terminal T IN with cascade connection, n multipliers 3-1 
to 3-n for multiplying signal input to the input terminal 
T IN and output signals of respective delay units 2-1 to 




2-(n-l) by filter coefficients h (0) to h <n-l) 
respectively, and an adder 4 for adding n output signals 
of the multipliers 3-1 to 3-n to output to an output 
terminal T OUT . 

As typical method for designing such the linear- 
phase FIR filter , for instance , these is known Remez 
Exchange algorithms which Parks, T.W. and McCleLLan, J.H. 
et al. apply it to the linear-phase FIR filter (Parks, 
T.W. and McClellan, J.H.: "Chebyshev Approximation for 
Non-recursive Digital Filters with Linear Phase", IEEE 
Trans. Circuit Theory, CT-19, 2 ,pp . 189-194 , 1972, as well 
as Rabiner, L.R., McClellan, J.H. and Parks, T.W.: "FIR 
Digital Filter Design Techniques Using Weighted Chebyshev 
Approximation", Proc . IEEE, Vol 63, April, pp. 595-610, 
1975) . 

The Remez Exchange algorithms is an algorithms in 
which a weighted approximation error is approximated such 
that the weighted approximation error is made to 
configure equi-ripple to the desired amplitude 
characteristics . 

There is known a resolution conversion of picture 
that utilizes the sampling rate conversion as the 
application of the filter processing using the linear- 
phase FIR filter. 

In this resolution conversion, multi-rate filter 



which has an interpolator, a decimeter and the linear- 
phase FIR filter as element- technique is employed (see, 
P.P.Vaidyanathan: "Multirate System and Filter Banks", 
Prentice Hall, 1992) . 
5 In use of the multi-rate filter, generally, the 

linear-phase FIR filter is made to use in such a way as 
to execute polyphase-sort (dissolution) in order to 
adjust the interpolator. Both the interpolator and the 
decimeter form periodic time invariance systems, thus 
10 having different characteristics from the time invariance 
system . 

Distortion so called as the chessboard distortion 
on lattice occurs in the resolution conversion of the 
picture caused by the periodic time invariance property 
15 of the interpolator . 

Accordingly, Harada , and Takaie considered the 
condition for avoiding such chessboard distortion from 
the viewpoint of zero-point arrangement of filter (see, 
Yasuhiro Harada, Hitoshi Kiya, : "Multi-rate Filter 
20 without Accompanying Chessboard Distortion and its Zero- 
point Arrangement" The technical Report of IEICE CAS96- 
78, ppl-6, 1997-01) . 

Transfer function H (z) for the multi-rate filter 
without accompanying the chessboard distortion will be 
25 discussed. The transfer function H (z) is capable of 



being found in such a way as to multiply a transfer 
function K (z) of the linear-phase FIR filter 
(hereinafter referred to as an equalizer) designed by a 
method in some kind by the transfer function Z (z) of the 
5 zero-point in order to avoid the chessboard distortion 
later . 



H(z)= Z{z)" K{z) (l) 

10 zH=Uz-' + r ! f + z- M <2) 

Here, a linear-phase FIR filter fixed beforehand 
such as the transfer function Z (z) of the zero-point for 
avoiding the chessboard distortion is called as a pre- 
15 filter. 

Fig. 2 shows an example of a frequency response of 
the multi-rate filter and weighted approximation error, 
in which the chessboard distortion is avoided by 
multiplying the equalizer designed by the use of Remez 
20 Exchange algorithms and the pre-filter. 

However, it. suffers from the disadvantage in the 
avoiding method of the chessboard distortion according to 
the aforementioned method. 

Namely, as illustrated in Fig. 2C, in the multi- 
25 rate filter having the transfer function H (z) designed 
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depending on the conventional method, the equi-ripple of 
the weighted approximation error designed depending on 
the Remez Exchange algorithms is not established. 

Further, as illustrated in Fig. 2B, the multi- 
5 rate filter designed depending on the conventional 

method, has gain of pass band with non-fixed value in 
which right end is attenuated. 

If the resolution conversion is executed using 
such filter, contours of picture appears fuzzy, thus is 
10 adversely affected on equality of the picture. 

The attenuation of this pass band can not be 
avoided even though the number of coefficient of filter 
is increased. 

15 SUMMARY OF THE INVENTION 

An object of the present invention is to* provide an 
FIR filter and setting method of its coefficients in 
which the equi-ripple of the weighted approximation error 
is not established, and further it is possible to 

20 maintain a gain of the pass band at approximately 
constant value. 

In order to achieve the above-mentioned objects, an 
FIR filter of the present invention has a configuration 
in which an impulse response is expressed by using a 

25 finite time length, this impulse response is equivalent 



to a filter coefficient of the FIR filter, and the FIR 
filter' s transfer function H (z) is related to a transfer 
function Z (z) of a pre-filter, wherein the filter 
coefficient is set by performing a weighted approximation 
5 to the desired characteristics in relation to the 
frequency response of the pre-filter. 

Also, an FIR filter of the present invention has a 
configuration in which the FIR filter's impulse response 
is expressed by using a finite time length, this impulse 
10 response is equivalent to a filter coefficient of the FIR 
filter, and whose transfer function H (z) is related to a 
transfer function Z (z) of a pre-filter, and a transfer 
function K (z) of an equalizer, wherein the filter 
coefficient is set on the basis of an amplitude 
15 characteristic of the equalizer that is obtained in such 
a way as to execute weighted approximation to the desired 
characteristics in relation to a frequency response of 
the pre-filter. 

Also, a setting method of filter coefficients of an 
20 FIR filter of the present invention has a method in which 
the FIR filter' s impulse response is expressed by using a 
finite time length, this impulse response is equivalent 
to a filter coefficient of the FIR filter, and whose 
transfer function H (z) is related to a transfer function 
25 Z (z) of a pre-filter, wherein the filter coefficient is 
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calculated by performing a weighted approximation to the 
desired characteristics in relation to a frequency 
response of the pre-filter. 

Also, a setting method of filter coefficients of an 
5 FIR filter of the present invention has a method in which 
the FIR filter's impulse response is expressed by using a 
finite time length, this impulse response is equivalent 
to a filter coefficient of the FIR filter, and whose 
transfer function H (z) is related to a transfer function 

10 Z (z) of a pre-filter and a transfer function K (z) of an 
equalizer, wherein the filter coefficient is calculated 
depending on an amplitude characteristic of the 
equalizer, which is obtained in such a way as to execute 
weighted approximation to the desired characteristics in 

15 relation to a frequency response of the pre-filter. 

Also, a setting method of filter coefficients of an 
FIR filter of the present invention has a method in which 
weighted approximation is executed to the desired 
characteristics by using the Remez Exchange algorithms 

20 taking into account a frequency response of a pre-filter. 

Also, a setting method of filter coefficients of an 
FIR filter of the present invention has a method in which 
the FIR filter's impulse response is expressed by a 
finite time length, and this impulse response is 

25 equivalent to a filter coefficient, the setting method of 
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a filter coefficient of the FIR filter comprises: a. first 
step for generating an interpolation polynomial equation 
for interpolating an amplitude characteristic from an 
extreme value point of the amplitude characteristic of a 
5 frequency; a second step for determining a new extreme 
value point from the amplitude characteristic obtained 
from the interpolation polynomial equation that is 
generated in the first step; a third step for judging 
whether or not a position of the extreme value is 
10 approximated within required range by repeating the 

operation in the first step and the second step; and a 
fourth step for finding the filter coefficient from the 
approximated amplitude characteristic obtained in the 
third step. 

15 Also, a setting method of filter coefficients of an 

FIR filter of the present invention has a method in which 
* there is further provided an initial setting step for 
carrying out, at least, setting of the FIR filter, 
setting of the band, setting of coefficient of a pre- 

20 filter, and setting of initial extreme value point, 
before executing the operation the first step. 

Also, a setting method of filter coefficients of an 
FIR filter of the present invention has a method, in the 
second step and the third step, the extreme-value of 

25 weighted approximation error calculated from the extreme- 
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value point used for the interpolation is searched for 
the entire approximation band, then obtained extreme- 
value is taken to be a new extreme- value point, and it is 
judged that the optimum approximation is obtained when 
5 the position of the extreme-value is not changed. 

Also, a setting method of filter coefficients of an 
FIR filter of the present invention has a method, in the 
fourth step, the filter coefficient is calculated by 
performing a weighted approximation to the desired 
10 characteristics in relation to a frequency response of 
the pre-filter. 

Also, a setting method of filter coefficients of an 
FIR filter of the present invention has a method, in the 
fourth step, the filter coefficient is calculated 
15 depending on amplitude characteristic of the equalizer 
obtained in such a way as to execute weighted 
approximation to the desired characteristics in relation 
to a frequency response of the pre-filter. 

Also, an FIR filter whose impulse response is 
20 expressed by using a finite time length, the impulse 

response being equivalent to a filter coefficient of the 
FIR filter, the FIR filter having arbitrary number of 
taps, and whose transfer function H (z) is related to a 
transfer function Z (z) of a pre-filter, wherein the 
25 filter coefficient is set in such a way as to perform a 
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weighted approximation to the desired characteristics so 
as to satisfy attenuation quantity of a stop band in 
relation to a frequency response of the pre-filter which 
satisfies the attenuation quantity of the stop band, when 
5 the number of taps is variable and the bandwidth is 
fixed . 

Also, an FIR filter whose impulse response is 
expressed by using a finite time length, the impulse 
response being equivalent to a filter coefficient of the 

10 FIR filter, the FIR filter having arbitrary number of 

taps, and whose transfer function H (z) is related to a 
transfer function Z (z) of a pre-filter and transfer 
function K (z) of an equalizer, wherein the filter 
coefficient is set on the basis of an amplitude 

15 characteristic of the equalizer obtained in such a way as 
to perform a weighted approximation to the desired 
characteristics so as to satisfy an attenuation quantity 
of a stop band in relation to frequency response of the 
pre-filter which satisfies the attenuation quantity of 

20 the stop band, when the number of taps is variable and 
the bandwidth is fixed. 

Also, a setting method of filter coefficients of an 
FIR filter of the present invention has a method in which 
the FIR filter's impulse response is expressed by using a 

25 finite time length, this impulse response is equivalent 
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to a filter coefficient of the FIR filter, the FIR filter 
has arbitrary number of tap, and whose transfer function 
H (z) is related to a transfer function Z (z) of a pre- 
filter, wherein the filter coefficient is calculated by 
5 performing a weighted approximation to the desired 

characteristics so as to satisfy attenuation quantity of 
a stop band in relation to the frequency response of the 
pre-filter which satisfies attenuation quantity of the 
stop band, when band is made to fix and the number of tap 
10 is made variable . 

Also, a setting method of filter coefficients of an 
FIR filter of the present invention has a method in which 
the FIR filter's impulse response is expressed by using a 
finite time length, this impulse response is equivalent 
15 to a filter coefficient of the FIR filter, the FIR filter 
has arbitrary number of tap, and whose transfer function 
H (z) is related to a transfer function Z (z) of a pre- 
filter and a transfer function K (z) of an equalizer, 
wherein the filter coefficient is calculated depending on 
20 an amplitude characteristic of the equalizer obtained in 
such a way as to execute weighted approximation to the 
desired characteristics so as to satisfy attenuation 
quantity of a stop band in relation to the frequency 
response of the pre-filter which satisfies attenuation 
25 quantity of the stop band, when band is made to fix and 



the number of tap is made variable. 

Also, a, setting method of filter coefficients of an 
FIR filter of the present invention has a method in which 
the FIR filter's impulse response is expressed by using a 
5 finite time length, this impulse response is equivalent 
to a filter coefficient of the FIR filter, and the FIR 
filter's number of tap is variable, and whose band is 
fixed, the setting method of a filter coefficient of an 
FIR filter comprises the steps of: a first step for 

10 generating interpolation polynomial equation for 

interpolating an amplitude characteristic from an extreme 
value point of the amplitude characteristic of the 
frequency; a second step for determining new extreme 
value point from the amplitude characteristic obtained 

15 from the interpolation polynomial equation that is 

obtained in the first step; a third step for judging 
whether or not position of the extreme value is 
approximated within required range while repeating the 
first step and the second step; a fourth step for 

20 examining attenuation quantity of a stop band from the 
approximated amplitude characteristic obtained in the 
third step; a fifth step for comparing the examined 
attenuation quantity with the attenuation quantity of the 
stop band thus designated to judge whether or not result 
25 of the comparison satisfies predetermined condition; a 
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sixth step for changing the number of tap when result of 
the comparison of the fifth step does not satisfy the 
predetermined condition; and a seventh step for finding 
the filter coefficient from the amplitude characteristic 
5 thus approximated depending on the third step which 

satisfies the predetermined condition in the fifth step. 

Also a setting method of filter coefficients of an 
FIR filter of the present invention has a method , wherein 
there is provided at least an initial setting step for 
10 carrying out setting of the FIR filter , setting of the 

band, setting of coefficient of a pre-filter, and setting 
of initial extreme value point, before executing the 
first step. 

Also, a setting method of filter coefficients of an 
15 FIR filter of the present invention has a method, in the 
fourth step, the minimum attenuation quantity in the stop 
band is examined, and in the sixth step, the number of 
the tap is increased. 

Also, a setting method of filter coefficients of an 
20 FIR filter of the present invention has a method, in the 
above described seventh step, the filter coefficient is 
calculated by performing a weighted approximation with 
reference to the desired characteristics so as to satisfy 
attenuation quantity of a stop band in relation to 
25 frequency response of the pre-filter that satisfies the 
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attenuation quantity of the stop band when a band is 
fixed and the number of tap is made variable. 

Also, a setting method of filter coefficients of an 
FIR filter of the present invention has a method, in the 
5 above described seventh step, the filter coefficient is 
calculated depending on an amplitude characteristic of an 
equalizer obtained in such a way as to execute weighted 
approximation with reference to the desired 

characteristics so as to satisfy attenuation quantity of 

10 a stop band in relation to frequency response of the pre- 
f ilter that satisfies the attenuation quantity of the 
stop band when a band is fixed and the number of tap is 
made variable . 

Also, an FIR filter of the present invention has a 

15 configuration in which the FIR filter' s impulse response 
is expressed by using a finite time length, this impulse 
response is equivalent to a filter coefficient of the FIR 
filter, the FIR filter has arbitrary number of tap, and 
whose transfer function H (z) is related to a transfer 

20 function Z (z) of a pre-f ilter, wherein the filter 

coefficient is set in such a way as to execute weighted 
approximation to the desired characteristics so as to 
satisfy attenuation quantity of a stop band in relation 
to frequency response of the pre-filter which satisfies 

25 attenuation quantity of a stop band, when the number of 
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taps is made to fix and band setting is changeable. 

Also, an FIR filter of the present invention has a 
configuration in which the FIR filter' s impulse response 
is expressed by using a finite time length, this impulse 
5 response is equivalent to a filter coefficient of the FIR 
filter, the FIR filter has arbitrary number of tap, and 
whose transfer function H (z) is related to a transfer 
function Z (z) of a pre-filter and transfer function K 
(z) of an equalizer, wherein the filter coefficient is 
10 set on the basis of an amplitude characteristic of the 
equalizer obtained in such a way as to execute weighted 
approximation to the desired characteristics so as to 
satisfy attenuation quantity of a stop band in relation 
to frequency response of the pre-filter which satisfies 
15 attenuation quantity of a stop band, when the number of 
taps is made to fix and band setting is changeable. 

Also, a setting method of filter coefficients of an 
FIR filter of the present invention has a method in which 
the FIR filter' s impulse response is expressed by using a 
20 finite time length, this impulse response is equivalent 

to a filter coefficient of the FIR filter, the FIR filter 
has arbitrary number of tap, and whose transfer function 
H (z) is related to a transfer function Z (z) of a pre- 
filter, wherein the filter coefficient is calculated by 
25 performing a weighted approximation to the desired 
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characteristics so as to satisfy attenuation quantity of 
a stop band in relation to the frequency response of the 
pre-filter which satisfies attenuation quantity of the 
stop band, when the number of taps is made to fix and 
5 band setting is changeable. 

Also, a setting method of filter coefficients of an 
FIR filter of the present invention has a method in which 
the FIR filter' s impulse response is expressed by using a 
;g finite time length, this impulse response is equivalent 

J 10 to a filter coefficient of the FIR filter, the FIR filter 

!H has arbitrary number of tap, and whose transfer function 

|_ H (z) is related to a transfer function Z (z) of a pre- 

£J filter and a transfer function K (z) of an equalizer, 

tl wherein the filter coefficient is calculated depending on 

15 an amplitude characteristic of the equalizer obtained in 
such a way as to execute weighted approximation to the 
desired characteristics so as to satisfy attenuation 
quantity of a stop band in relation to the frequency 
response of the pre-filter which satisfies attenuation 
20 quantity of the stop band, when the number of taps is 
made to fix and band setting is changeable . 

Also, a setting method of filter coefficients of an 
FIR filter of the present invention has a method in which 
the FIR filter' s impulse response is expressed by using a 
25 finite time length, this impulse response is equivalent 
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to a filter coefficient of the FIR filter, and the FIR 
filter' s number of tap is fixed, and whose band setting 
is changeable, the setting method of a filter coefficient 
of an FIR filter comprises the steps of: a first step for 
5 generating interpolation polynomial equation for 

interpolating an amplitude characteristic from an extreme 
value point of the amplitude characteristic of the 
frequency; a second step for determining new extreme 
value point from the amplitude characteristic obtained 
10 from the interpolation polynomial equation that is 

obtained in the first step; a third step for judging 
whether or not position of the extreme value is 
approximated within required range while repeating the 
first step and the second step; a fourth step for 
15 examining attenuation quantity of a stop band from the 
approximated amplitude characteristic obtained in the 
third step; a fifth step for comparing the examined 
attenuation quantity with the attenuation quantity of the 
stop band thus designated to judge whether or not result 
20 of the comparison satisfies predetermined condition; a 
sixth step for changing the band setting when result of 
the comparison of the fifth step does not satisfy the 
predetermined condition; and a seventh step for finding 
the filter coefficient from the amplitude characteristic 
25 thus approximated depending on the third step which 
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satisfies the predetermined condition in the fifth step. 

Also, a setting method of filter coefficients of an 
FIR filter of the present invention has a method, wherein 
there is provided an initial setting step for carrying 
5 out, at least, setting of the FIR filter, setting of the 
band, setting of coefficient of a pre-filter, and setting 
of initial extreme value point, before executing the 
;3 first step. 

l5 Also, a setting method of filter coefficients of an 

:; p 10 FIR filter of the present invention has a method, in the 

ffl fourth step, the minimum attenuation quantity in the stop 

band is examined. 
;:L? Also, a setting method of filter coefficients of an 

^ FIR filter of the present invention, in the above 

15 described seventh step, the filter coefficient is 

calculated by performing a weighted approximation with 
reference to the desired characteristics so as to satisfy 
attenuation quantity of a stop band in relation to 
frequency response of the pre-filter that satisfies the 
20 attenuation quantity of the stop band when the number of 
taps is fixed and the band setting is changeable . 

Also, a setting method of filter coefficients of an 
FIR filter of the present invention has a method, in the 
above described seventh step, the filter coefficient is 
25 calculated depending on an amplitude characteristic of an 
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equalizer obtained in such a way as to execute weighted 
approximation with reference to the desired 
characteristics so as to satisfy attenuation quantity of 
a stop band in relation to frequency response of the pre- 
5 filter that satisfies the attenuation quantity of the 
stop band when the number of taps is fixed and the band 
setting is changeable. 

Also, an FIR filter of the present invention has a 
configuration in which the FIR filter's impulse response 

10 is expressed by using a finite time length, this impulse 
response is equivalent to a filter coefficient of the FIR 
filter, the FIR filter has arbitrary number of tap, and 
whose transfer function H (z) is related to a transfer 
function Z (z) of a pre-filter, wherein the filter 

15 coefficient is set in such a way as to execute weighted 
approximation to the desired characteristics so as to 
satisfy attenuation quantity of a stop band in relation 
to frequency response of the pre-filter which satisfies 
attenuation quantity of a stop band, when the number of 

20 taps is variable and band setting is changeable. 

Also, an FIR filter of the present invention has a 
configuration in which the FIR filter's impulse response 
is expressed by using a finite time length, this impulse 
response is equivalent to a filter coefficient of the FIR 

25 filter, the FIR filter has arbitrary number of tap, and 
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whose transfer function H (z) is related to a transfer 
function Z (z) of a pre-filter and transfer function K 
(z) of an equalizer, wherein the filter coefficient is 
set on the basis of an amplitude characteristic of the 
5 equalizer obtained in such a way as to execute weighted 
approximation to the desired characteristics so as to 
satisfy attenuation quantity of a stop band in relation 
to frequency response of the pre-filter which satisfies 
attenuation quantity of a stop band, when the number of 

10 taps is variable and band setting is changeable. 

Also, a setting method of filter coefficients of an 
FIR filter of the present invention has a way in which 
the FIR filter's impulse response is expressed by using a 
finite time length , this impulse response is equivalent 

15 to a filter coefficient pf the FIR filter, the FIR filter 
has arbitrary number of tap, and whose transfer function 
H (z) is related to a transfer function Z (z) of a pre- 
filter, wherein the filter coefficient is calculated by 
performing a weighted approximation to the desired 

20 characteristics so as to satisfy attenuation quantity of 
a stop band in relation to the frequency response of the 
pre-filter which satisfies attenuation quantity of the 
stop band, when the number of taps is variable and band 
setting is changeable. 

25 Also, a setting method of filter coefficients of an 
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FIR filter of the present invention has a way in which 
the FIR filter' s impulse response is expressed by using a 
finite time length, this impulse response is equivalent 
to a filter coefficient of the FIR filter, the FIR filter 
5 has arbitrary number of tap, and whose transfer function 
H (z) is related to a transfer function Z (z) of a pre- 
filter and a transfer function K (z) of an equalizer, 
wherein the filter coefficient is calculated depending on 
an amplitude characteristic of the equalizer obtained in 

10 such a way as to execute weighted approximation to the 
desired characteristics so as to satisfy attenuation 
quantity of a stop band in relation to the frequency 
response of the pre-filter which satisfies attenuation 
quantity of the stop band, when the number of taps is 

15 variable and band setting is changeable. 

Also, a setting method of filter coefficients of an 
FIR filter of the present invention has a way in which 
the FIR filter 7 s impulse response is expressed by using a 
finite time length, this impulse response is equivalent 

20 to a filter coefficient of the FIR filter, and the FIR 
filter's number of tap is variable, and whose band 
setting is changeable, the setting method of a filter 
coefficient of an FIR filter comprises the steps of: a 
first step for generating interpolation polynomial 

25 equation for interpolating an amplitude characteristic 
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from an extreme value point of the amplitude 
characteristic of the frequency; a second step for 
determining new extreme value point from the amplitude 
characteristic obtained from the interpolation polynomial 
5 equation that is obtained in the first step; a third step 
for judging whether or not position of the extreme value 
is approximated within required range while repeating the 
first step and the second step; a fourth step for 
examining attenuation quantity of a stop band from the 

10 approximated amplitude characteristic obtained in the 
third step; a fifth step for comparing the examined 
attenuation quantity with the attenuation quantity of the 
stop band thus designated to judge whether or not result 
of the comparison satisfies predetermined condition; a 

15 sixth step for changing the band setting when result of 
the comparison of the fifth step does not satisfy the 
predetermined condition; a seventh step for judging 
whether or not the present number of taps are capable of 
satisfying the attenuation quantity of the stop band 

20 after changing of the band in the sixth step; an eighth 
step for changing the number of taps when judgement is 
performed that the present number of taps do not satisfy 
the attenuation quantity of the stop band in the seventh 
step; and a ninth step for finding the filter coefficient 

25 from the amplitude characteristic thus approximated 
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depending on the third step which satisfies the 
predetermined condition in the fifth step. 

Also, a setting method of filter coefficients of an 
FIR filter of the present invention has a way, wherein 
5 there is provided at least an initial setting step for 
carrying out setting of the FIR filter , setting of the 
band, setting of coefficient of a pre-filter, and setting 
of initial extreme value point, before executing the 
=3 first step. 

if 10 Also, a setting method of filter coefficients of an 

=!S FIR filter of the present invention has a way, in the 

fourth step, the minimum attenuation quantity in the stop 
i;n band is examined, and in the eighth step, the number of 

i!f§ the tap is increased. 

M 15 Also, a setting method of filter coefficients of an 

FIR filter of the present invention has a way, in the 
above described ninth step, the filter coefficient is 
calculated by performing a weighted approximation with 
reference to the desired characteristics so as to satisfy 

20 attenuation quantity of a stop band in relation to 

frequency response of the pre-filter that satisfies the 
attenuation quantity of the stop band when the number of 
taps is variable and band setting is changeable. 

Also, a setting method of filter coefficients of an 

25 FIR filter of the present invention has a way, in the 
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above described ninth step, the filter coefficient is 
calculated depending on an amplitude characteristic of an 
equalizer obtained in such a way as to execute weighted 
approximation with reference to the desired 
5 characteristics so as to satisfy attenuation quantity of 
a stop band in relation to frequency response of the pre- 
filter that satisfies the attenuation quantity of the 
stop band when the number of taps is variable and the 
band setting is changeable. 

10 Also, an FIR filter of the present invention has a 

configuration in which the FIR filter' s impulse response 
is expressed by using a finite time length, this impulse 
response is equivalent to a filter coefficient of the FIR 
filter , the FIR filter has arbitrary number of tap, and 

15 whose transfer function H (z) is related to a transfer 
function Z (z) of a pre-f ilter , wherein the filter 
coefficient is set in such a way as to execute weighted 
approximation to the desired characteristics so as to 
satisfy attenuation quantity of a stop band in relation 

20 to frequency response of the pre-f ilter. through which the 
attenuation quantity of the designated frequency of a 
transition band is passed, and which satisfies 
attenuation quantity of a stop band, when the number of 
taps is made to fix and band setting is changeable. 

25 Also, an FIR filter of the present invention has a 
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configuration in which the FIR filter' s impulse response 
is expressed by using a finite time length, this impulse 
response is equivalent to a filter coefficient of the FIR 
filter, the FIR filter has arbitrary number of tap, and 
5 whose transfer function H (z) is related to a transfer 
function Z (z) of a pre-filter and transfer function K 
(z) of an equalizer, wherein the filter coefficient is 
set on the basis of an amplitude characteristic of the 
equalizer obtained in such a way as to execute weighted 

10 approximation to the desired characteristics so as to 

satisfy attenuation quantity of a stop band in relation 
to frequency response of the pre-filter through which the 
attenuation quantity of the designated frequency of a 
transition band is passed, and which satisfies 

15 attenuation quantity of a stop band, when the number of 
taps is made to fix and band setting is changeable. 

Also, a setting method of filter coefficients of an 
FIR filter of the present invention has a way in which 
the FIR filter' s impulse response is expressed by using a 

20 finite time length, this impulse response is equivalent 

to a filter coefficient of the FIR filter, the FIR filter 
has arbitrary number of tap, and whose transfer function 
H (z) is related to a transfer function Z (z) of a pre- 
filter, wherein the filter coefficient is calculated by 

25 performing a weighted approximation to the desired 
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characteristics so as to satisfy attenuation quantity of 
a stop band in relation to the frequency response of the 
pre-filter through which attenuation quantity of the 
designated frequency of a stop band is passed, and which 
5 satisfies attenuation quantity of the stop band, when the 
number of taps is made to fix and band setting is 
changeable . 

.«! Also, a setting method of filter coefficients of an 

FIR filter of the present invention has a way in which 
:: p 10 the FIR filter' s impulse response is expressed by using a 

in finite time length, this impulse response is equivalent 

;^ to a filter coefficient of the FIR filter, the FIR filter 

;L; has arbitrary number of tap, and whose transfer function 

;:L; H (z) is related to a transfer function Z (z) of a pre- 

15 filter and a transfer function K (z) of an equalizer, 

wherein the filter coefficient is calculated depending on 
an amplitude characteristic of the equalizer obtained in 
such a way as to execute weighted approximation to the 
desired characteristics so as to satisfy attenuation 
20 quantity of a stop band in relation to the frequency 
response of the pre-filter through which attenuation 
quantity of the designated frequency of a stop band is 
passed, and which satisfies attenuation quantity of the 
stop band, when the number of taps is made to fix and 
25 band setting is changeable. 
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Also, a setting method of filter coefficients of an 
FIR filter of the present invention has a way in which 
the FIR filter' s impulse response is expressed by using a 
finite time length, this impulse response is equivalent 
5 to a filter coefficient of the FIR filter, and the FIR 
filter' s number of tap is fixed, and whose band setting 
is changeable, the setting method of a filter coefficient 
of an FIR filter comprises the steps of: a first step for 
generating interpolation polynomial equation for 

10 interpolating an amplitude characteristic from an extreme 
value point of the amplitude characteristic of the 
frequency; a second step for determining new extreme 
value point from the amplitude characteristic obtained 
from the interpolation polynomial equation that is 

15 obtained in the first step; a third step for judging 
whether or not position of the extreme value is 
approximated within required range while repeating the 
first step and the second step; a fourth step for 
examining attenuation quantity of a stop band from the 

20 approximated amplitude characteristic obtained in the 
third step; a fifth step for comparing the examined 
attenuation quantity in the fourth step with the 
attenuation quantity of the stop band thus designated to 
judge whether or not result of the comparison satisfies 

25 predetermined condition; a sixth step for changing the 
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h>and setting when result of the comparison of the fifth 
step does not satisfy the predetermined condition; a 
seventh step for examining attenuation quantity of the 
designated frequency of the transition band which 
5 attenuation quantity satisfies predetermined condition in 
the fifth step; an eighth step for comparing the 
attenuation quantity of the designated frequency of the 
transition band that is examined in the seventh step with 
the attenuation quantity of the designated transition 

10 band, and for judging whether or not result of comparison 
satisfies predetermined condition; a ninth step for 
changing setting of the band when result of comparison of 
the seventh step does not satisfy the predetermined 
condition; and a tenth step for finding the filter 

15 coefficient from the amplitude characteristic thus 
approximated depending on the seventh step which 
amplitude characteristic satisfies the predetermined 
condition. 

Also, a setting method of filter coefficients of an 
20 FIR filter of the present invention has a way, wherein 
there is provided at least an initial setting step for 
carrying out setting of the FIR filter, setting of the 
band, setting of coefficient of a pre-filter, and setting 
of initial extreme value point, before executing the 
25 first step. 
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Also, a setting method of filter coefficients of an 
FIR filter of the present invention has a way, in the 
fourth step, the minimum attenuation quantity in the stop 
band is examined. 
5 Also, a setting method of filter coefficients of an 

FIR filter of the present invention has a way, in the 
above described tenth step, the filter coefficient is 
calculated by performing a weighted approximation with 
reference to the desired characteristics so as to satisfy 
10 attenuation quantity of a stop band in relation to 

frequency response of the pre-filter that satisfies the 
attenuation quantity of the stop band, and that causes 
the attenuation quantity of the designated frequency of 
the transition band to pass when the number of taps is 
15 fixed and band setting is changeable. 

Also, a setting method of filter coefficients of an 
FJR filter of the present invention has a way, in the 
above described tenth step, the filter coefficient is 
calculated depending on an amplitude characteristic of an 
20 equalizer obtained in such a way as to execute weighted 
approximation with reference to the desired 

characteristics so as to satisfy attenuation quantity of 
a stop band in relation to frequency response of the pre- 
filter that satisfies the attenuation quantity of the 
25 stop band, and that causes the attenuation quantity of 
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the designated frequency of the transition band to pass 
when the number of taps is fixed and the band setting is 
changeable . 

Also, an FIR filter of the present invention has a 
5 configuration in which the FIR filter's impulse response 
is expressed by using a finite time length, this impulse 
response is equivalent to a filter coefficient of the FIR 
filter, the FIR filter has arbitrary number of tap, and 
whose transfer function H (z) is related to a transfer 

10 function Z (z) of a pre-filter, wherein the filter 

coefficient is set in such a way as to execute weighted 
approximation to the desired characteristics so as to 
satisfy attenuation quantity of a stop band in relation 
to frequency response of the pre-filter through which the 

15 attenuation quantity of the designated frequency of a 
transition band is passed, and which satisfies 
attenuation quantity of a stop band, when the number of 
taps is variable and band setting is changeable. 

Also, an FIR filter of the present invention has a 

20 configuration in which the FIR filter's impulse response 
is expressed by using a finite time length, this impulse 
response is equivalent to a filter coefficient of the FIR 
filter, the FIR filter has arbitrary number of tap, and 
whose transfer function H (z) is related to a transfer 

25 function Z (z) of a pre-filter and transfer function K 
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(z) of an equalizer, wherein the filter coefficient is 
set on the basis of an amplitude characteristic of the 
equalizer obtained in such a way as to execute weighted 
approximation to the desired characteristics so as to 
5 satisfy attenuation quantity of a stop band in relation 
to frequency response of the pre-filter through which the 
attenuation quantity of the designated frequency of a 
transition band is passed, and which satisfies 
attenuation quantity of a stop band, when the number of 

10 taps is variable and band setting is changeable. 

Also, a setting method of filter coefficients of an 
FIR filter of the present invention has a way in which 
the FIR filter' s impulse response is expressed by using a 
finite time length, this impulse response is equivalent 

15 to a filter coefficient of the FIR filter, the FIR filter 
has arbitrary number of tap, and whose transfer function 
H (z) is related to a transfer function Z (z) of a pre- 
filter, wherein the filter coefficient is calculated by 
performing a weighted approximation to the desired 

20 characteristics so as to satisfy attenuation quantity of 
a stop band in relation to the frequency response of the 
pre-filter through which attenuation quantity of the 
designated frequency of a stop band is passed, and which 
satisfies attenuation quantity of the stop band, when the 

25 number of taps variable and band setting is changeable. 
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Also, a setting method of filter coefficients of an 
FIR filter of the present invention has a way in which 
the FIR filter's impulse response is expressed by using a 
finite time length, this impulse response is equivalent 
5 to a filter coefficient of the FIR filter, the FIR filter 
has arbitrary number of tap, and whose transfer function 
H (z) is related to a transfer function Z (z) of a pre- 
f liter and a transfer function K (z) of an equalizer, 
wherein the filter coefficient is calculated depending on 

10 an amplitude characteristic of the equalizer obtained in 
such a way as to execute weighted approximation to the 
desired characteristics so as to satisfy attenuation 
quantity of a stop band in relation to the frequency 
response of the pre-filter through which attenuation 

15 quantity of the designated frequency of a stop band is 
passed, and which satisfies attenuation quantity of the 
stop band, when the number of taps is variable and band 
setting is changeable . 

Also, a setting method of filter coefficients of an 

20 FIR filter of the present invention has a way in which 

the FIR filter' s impulse response is expressed by using a 
finite time length, this impulse response is equivalent 
to a filter coefficient of the FIR filter, and the FIR 
filter's number of tap is variable, and whose band 

25 setting is changeable, the setting method of a filter 
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coefficient of an FIR filter comprises the steps of: a 
first step for generating interpolation polynomial 
equation for interpolating an amplitude characteristic 
from an extreme value point of the amplitude 
5 characteristic of the frequency; a second step for 

determining new extreme value point from the amplitude 
characteristic obtained from the interpolation polynomial 
equation that is obtained in the first step; a third step 
for judging whether or not position of the extreme value 
10 is approximated within required range while repeating the 
first step and the second step; a fourth step for 
examining attenuation quantity of a stop band from the 
approximated amplitude characteristic obtained in the 
third step; a fifth step for comparing the examined 
15 attenuation quantity in the fourth step with the 

attenuation quantity of the stop band thus designated to 
judge whether or not result of the comparison satisfies 
predetermined condition; a sixth step for changing the 
band setting when result of the comparison of the fifth 
20 step does not satisfy the predetermined condition; a 
seventh step for judging whether or not the present 
number of taps is capable of satisfying attenuation 
quantity of a stop band after changing of the band in the 
sixth step; an eighth step for changing the number of 
25 taps when judgement is performed that the present number 
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of taps can not satisfy the attenuation quantity in the 
seventh step; a ninth step for examining attenuation 
quantity of the designated frequency of the transition 
band which attenuation quantity satisfies predetermined 
5 condition in the fifth step; a tenth step for comparing 
the attenuation quantity of the designated frequency of 
the transition band that is examined in the ninth step 
with the attenuation quantity of the designated 
transition band, and for judging whether or not result of 
10 comparison satisfies predetermined condition; an eleventh 
step for changing setting of the band when result of 
comparison of the tenth step does not satisfy the 
predetermined condition; a twelfth step for judging 
whether or not the present number of taps causes the 
15 signal to pass the designated frequency of the stop band 
after changing the band in the eleventh step; a 
thirteenth step changing the number of taps when 
judgement is performed that the present number of taps 
does not enable the designated frequency to be passed in 
20 the twelfth step; and a fourteenth step for finding the 

filter coefficient from the amplitude characteristic thus 
approximated depending on the tenth step which amplitude 
characteristic satisfies the predetermined condition. 

Also, a setting method of filter coefficients of an 
25 FIR filter of the present invention has a way, wherein 
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there is provided at least an initial setting step for 
carrying out setting of the FIR filter, setting of the 
band, setting of coefficient of a pre-filter, and setting 
of initial extreme value point, before executing the 
5 first step. 

Also, a setting method of filter coefficients of an 
FIR filter of the present invention has a way, in the 
fourth step, the minimum attenuation quantity in the stop 
band is examined, and both in the eighth step and in the 

10 thirteenth step, the number of the tap is increased. 

Also, a setting method of filter coefficients of an 
FIR filter of the present invention has a way, in the 
above described fourteenth step, the filter coefficient 
is calculated by performing a weighted approximation with 

15 reference to the desired characteristics so as to satisfy 
attenuation quantity of a stop band in relation to 
frequency response of the pre-filter that satisfies the 
attenuation quantity of the stop band, and that causes 
the attenuation quantity of the designated frequency of 

20 the transition band to pass when the number of taps is 
variable and band setting is changeable. 

Also, a setting method of filter coefficients of an 
FIR filter of the present invention has a way, in the 
above described fourteenth step, the filter coefficient 

25 is calculated depending on an amplitude characteristic of 
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an equalizer obtained in such a way as to execute 
weighted approximation with reference to the desired 
characteristics so as to satisfy attenuation quantity of 
a stop band in relation to frequency response of the pre- 
5 filter that satisfies the attenuation quantity of the 
stop band, and that causes the attenuation quantity of 
the designated frequency of the transition band to pass 
when the number of taps is variable and the band setting 
is changeable . 

10 According to the present invention, for instance, 

depending on the initial setting, the setting of the 
linear-phase FIR filter, the setting of the band, the 
setting of the coefficient of the pre-filter, and the 
setting of the extreme- value point are carried out. 

15 Subsequently, an interpolation polynomial equation is 

generated for interpolating amplitude characteristic from 
the present extreme value point of the amplitude 
characteristic of the frequency. 

Next, a new extreme value point is determined from 

20 the amplitude characteristic obtained from the 

interpolation polynomial equation that is obtained in the 
first step. 

Judgement is performed whether or not position of the 
extreme value is approximated within required range while 
25 repeating the first step and the second step. 
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Then, tinding of the filter coefficient is performed 
from the approximated amplitude characteristic obtained. 

Thus, the FIR filter whose coefficient is set, is 
capable of possessing the weighted approximation error 
5 with equi-ripple, further in which gain of the pass band 
is maintained constant value. 



10 invention will become clearer from the following 

description of the preferred embodiments with reference 
to the accompanying drawings , in which : 

Fig. 1 is a view illustrating a transversal type 
circuit configuration of an FIR filter; 

15 Figs . 2A to 2C are views illustrating one example of 

frequency response and weighted approximation error in 
order to avoid the chessboard distortion in the 
conventional method ; 

Figs . 3A to 3D are views illustrating impulse response 

20 of four cases where the FIR filter possesses linear- 
phase ; 

Fig. 4 is a view illustrating Q(e ju ) and R to four 
cases of the linear-phase FIR filter; 

Fig. 5 is a view illustrating example of weighted 
25 Chebyshev approximation; 



BRIEF DESCRIPTION OF THE DRAWINGS 



These and other objects and features of the present 
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Fig. 6 is a flowchart of the Remez Exchange algorithms 
while taking into account frequency response of the pre- 
f liter according to the present invention; 

Figs . 7A to 7C are views for explaining determining 
5 method of new extreme value of the weighted approximation 
error E(e j ") ; 

Fig. 8 is a view illustrating frequency response of 
low-pass filter which avoids the chessboard distortion 
designed according to the present invention; 
10 Figs . 9A and 9B are views for comparing the frequency 

response of the low-pass filter designed by the 
conventional technique with the present invention 
respectively ; 

Fig. 10 is a view illustrating weighted approximation 
15 error to the filter designed according to the present 
invention ; 

Fig. 11 is a view illustrating flowchart of algorithms 
for finding filter which satisfies attenuation quantity 
of the stop band; 
20 Fig. 12 is a view illustrating parameter of algorithms 

for finding filter that satisfies the attenuation 
quantity of the stop band and has the maximum end 
frequency of the pass band; 

Fig. 13 is a view illustrating initial frequency of 
25 the dichotomizing method in the algorithms for finding 
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filter having the maximum end frequency of the pass band 
which filter satisfies the attenuation quantity of the 
stop band; 

Figs.l4A to 14C are views illustrating change of band 
setting in loop of the first time; 

Figs. ISA and 15B are views illustrating change of 
band setting in loop after the second time and later ; 

Fig. 16 is a view illustrating frequency response of 
the filter having the maximum end frequency of pass band 
which filter satisfies attenuation quantity of the stop 
band ; 

Fig. 17 is a view illustrating parameter in the 
algorithms for finding filter having the minimum starting 
frequency of the stop band which filter satisfies 
attenuation quantity of the stop band; 

Fig. 18 is a view illustrating initial frequency of 
the dichotomizing method in the algorithms for finding 
filter having the minimum starting frequency of the pass 
band which filter satisfies the attenuation quantity of 
the stop band; 

Figs.l9A to 19C are views illustrating change of band 
setting in loop of the first time; 

Figs. 2 OA and 20B are views illustrating change of 
band setting in loop after the second time and later; 

Fig. 21 is a view illustrating frequency response of 
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the filter having the minimum starting frequency of stop 
band which filter satisfies attenuation quantity of the 
stop band; 

Fig. 22 is a view illustrating flowchart of the 
algorithms for finding filter with the minimum number of 
tap that satisfies attenuation quantity of the stop band; 

Fig. 23 is a view illustrating frequency response of 
the filter with the minimum number of tap that satisfies 
attenuation quantity of the stop band; 

Fig. 24 is a view illustrating a flowchart of the 
algorithms for finding filter that satisfies the 
attenuation quantity of the stop band and that enable the 
signal to be passed through the frequency point of the 
transition band; 

Fig. 25 is a view illustrating algorithms (1) for 
finding filter that satisfies the attenuation quantity of 
the stop band and that enable the signal to be passed 
through the frequency point of the transition band; 

Fig. 2 6 is a view illustrating initial frequency of 
the dichotomizing method in the algorithms for finding 
filter that satisfies the attenuation quantity of the 
stop band and that enables the signal to be passed 
through the frequency point of the transition band; 

Figs.27A to 29C^are views illustrating change of band 
setting in loop of the first time; 
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Figs.28A and 28B are views illustrating change of 
band setting in loop after the second time and later; 

Fig. 29 is a view illustrating frequency response (1) 
of the filter that satisfies the attenuation quantity of 
5 the stop band and that enables the signal to be passed 
through the frequency point of the transition band; 

Fig. 30 is a view illustrating algorithms (2) for 
finding filter that satisfies the attenuation quantity of 
the stop band and that enable the signal to be passed 
10 through the frequency point of the transition band; 

Fig. 31 is a view illustrating initial frequency of 
the dichotomizing method in the algorithms for finding 
filter that satisfies the attenuation quantity of the 
stop band and that enables the signal to be passed 
15 through the frequency point of the transition band; 

Figs . 32A to 32C are views illustrating change of band 
setting in loop of the first time; 

Figs.33A and 33B are views illustrating change of 
band setting in loop after the second time and later ; 
20 Fig. 34 is a view illustrating frequency response (2) 

of the filter that satisfies the attenuation quantity of 
the stop band and that enables the signal to be passed 
through the frequency point of the transition band; 
Fig. 35 is a view illustrating flowchart of the 
25 algorithms for finding filter with the minimum number of 
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tap that satisfies attenuation quantity of the stop band; 

Fig. 36 is a view illustrating frequency response of 
the filter with the minimum number of tap that satisfies 
attenuation quantity of the stop band; 

Fig. 37 is a view illustrating frequency response of 
the filter with the minimum number of tap that satisfies 
attenuation quantity of the stop band; 

Fig. 38 is a view illustrating flowchart of the 
algorithms for finding filter with the minimum number of 
tap that satisfies the attenuation quantity of the stop 
band and that enables the signal to be passed through the 
frequency band of the transition band; and 

Fig. 39 is a diagram illustrating frequency response 
of the filter with the minimum number of taps which 
filter satisfies the attenuation quantity of the stop 
band, and which filter enables the signal to be passed 
through the frequency point of the transition band. 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 

Below, preferred embodiments will be described with 

reference to the accompanying drawings . 

Preferred embodiments of the present invention will 

be explained with reference to the accompanying drawings 

below . 

A linear-phase FIR filter according to the present 
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invention is capable of equivalently adopting transversal 
type circuit configuration as illustrated in Fig . 1 . The 
FIR filter includes n-JL series-connected unit time delay 
elements 2-1 to 2-n-l, n multipliers 3-1 to 3-n having 
5 filter coefficients h(0) to h(n-l), where (n-1) 

multipliers 3-1 to 3-n-l are connected to input terminals 
of the corresponding unit time delay elements 2-1 to 2-n- 
1 and n-th multiplier 3-n is connected to an output 
terminal of n-th time unit time delay element 2-n-l, and 
10 an adder 4 connected to output terminals of the n 

multipliers 3-1 to 3-n. However, filter coefficients h, 
as described later in detail, can be obtained in such a 
way wherein the Remez Exchange algorithms is extended, 
and the desired amplitude characteristic is subjected to 
15 Chebyshev approximation under taking into account a 
frequency response of a pre-f liter , and the filter 
coefficients h are obtained from the approximated 
amplitude characteristics . 

Below, a concrete method for setting coefficients of 
20 the linear-phase FIR filter according to the present 

invention will be explained successively with reference 
to the drawings . 

In equation (3) , the transfer function H (z) of the 
linear-phase FIR filter having N-tap is defined as a 
25 filter which calculates the product of the transfer 
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function Z (z) of the pre-filter and the transfer 
function of the equalizer. 

H{z)^Z{zYK\z) -{3) 

5 

Here, the pre-filter is a U tap of the linear-phase 
FIR filter, and the equalizer is a N - (U - 1) tap 
linear-phase FIR filter respectively. The transfer 
function of the pre-filter is previously given. 

10 A design of a filter of the transfer function H (z) 

is the determination of the transfer function K (z) of 
the equalizer having N - (U - 1) taps such that the 
amplitude characteristic H (e ju ) is approached to the 
desired amplitude characteristic D (e jtJ ) . 

15 in the present embodiment, in order to design the 

transfer function H (z) under consideration of the 
amplitude characteristic K (e ju ) of such the pre-filter, 
a method in which the Remez Exchange algorithms is 
extended for solving Chebyshev approximation problem, is 

20 adopted. 

The number of taps allocated to the equalizer of the 
transfer function K (z) is defined as L = N - (U - 1) . 

The transfer function K (z) of the linear-phase FIR 
filter, as illustrated in Figs. 3A to 3D, is categorized 
25 with four cases because the transfer function thereof has 
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linear-phase. 

Concretely, it is categorized into four cases, case 1 
: odd number of taps, and even symmetry as illustrated in 



5 as illustrated in Fig.3B, case 3 : odd number of taps, 
and odd symmetry as illustrated in Fig.3C, and case 4 : 
even number of taps , and odd symmetry . 

With reference to the case 1 , the amplitude 
characteristic function K (e ju ) is made to leave as it 
10 is, and for the other cases 2 to 4 , rewrite as follows. 



Fig. 3A, case 2 



even number of taps , and even symmetry 



U-0/2 




Case 1 



(4-1) 



Case 2 




(4-2) 



U-l)/2 




15 



Case 3 



U-3)/2 



(4-3) 




Case 4 




(4-4) 
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Namely, the amplitude characteristic function K (e 3 ") 
is expressed by the product of the function Q (e 3 ") of 
the fixed parameters as illustrated in Fig. 4 and a cosine 
5 series P (e Dtt ) including the design parameters. 

Hereinafter, the upper limit of sum of the respective 
equations (4-1) to (4-4) is expressed as R - 1 - Namely, R 
is calculated as shown in Fig. 4. 

Further, a(n}, b(n}, c(n}, d(n) are expressed generally 

10 as P (n) . 

When the desired amplitude characteristic is taken to 
be D (e ju ) , and weight to respective frequencies are 
taken to be W (e jw ) , the weighted approximation error is 
defined as follows : 

15 

E{e jt " ) = w{e JW \de j0J ) - H{e jo) )} ( 5 ) 

H(e j(0 ) = I<{e ico ) • z{e jw ) = (far ) • l{e Jcu ) • z{e io> ) ( 6 ) 

20 Substitute equation (6) into equation (5) , the 

following equation is obtained. 



(7) 



47 



However 



, w(e J(0 \ fie 1 ") 



are expressed as follows 



Me*'" ) = u{e iM ) • de j< ° ) m z(e ico ) (8 ) 



w - \ de im ) ■ 



The equation (7) expresses the weighted approximation 
error of the linear-phase FIR filter in the four cases of 
case 1 to case 4 . 

The weighted Chebyshev approximation problem is to 

determine a(n\ b{n\ c(?t\ d(n) of the equations (4-1) to (4-4) 

that makes to bring the maximum value of | within 

designated frequency band into minimum, in the equation 
(5) . 

Hereinafter, there will be explained in relation to 
the concrete example . 

Here, the amplitude characteristic D (e ju ) is defined 
as illustrated in Fig. 5 and the following. 

= 1 (error within ±6 X , 0<u<u p ) 
C{e j '") = 0 (error within ±5 2 , co g <co<n) (10) 
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Note, when R is given, values of 5 1; 5 2 can not be 
designated arbitrarily, but it is possible to designate 

its ratio . 

The W (e^) is taken to be a constant value Wl in the 
pass band, and to be W2 in the stop band, thus being 
selected in such a way that Wl 5 X = W2 5 2 stands . For 
instance, Wl = 1, W2 = 5 X / 5 2 are made to select. At 
this time, next alternating theorem stands. 
Theorem 

The necessary and sufficient condition in order that 
the cosine series P (•*•> of <R - D order is the best 
weighted Chebyshev approximation to the target 
characteristic within the section (0, n) of co is as 
follows : 

(1) E (e iu ) obtains the extreme value at least (R + 
1) times within the section (0, n) . At this time, 
frequency within which the extreme value is obtained is 
taken to be . 0 < « 1< • 2 < ■ • < » <*-!> < u R " 

(2) The sign of the neighboring extreme values are 
different, and the absolute values of the whole extreme 
values are the same. Namely, the following conditions are 
satisfied. 



. , R - 1) 
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|4<? y, "]=|4 e;ai+ ') <i = °' i. r - x) <Ll) 

Consequently, )| is equal to the maximum value 

of Jf^e-*") within the section. 
5 There is known Remez Exchange algorithms on the basis 

of the alternating theorem as the technique for obtaining 
the best Chebyshev approximation (Rabiner , L.R- , 
McClellan, J.H. and Parks, T.W.: "FIR Digital Filter 
Design Techniques Using Weighted Chebyshev 
10 Approximation", Proc . IEEE, Vol 63, April, pp. 595-610, 
1975) . 

The Remez Exchange algorithms is an algorithms in 
which the desired amplitude characteristic is carried out 
Chebyshev approximation within frequency domain, and 
15 coefficients of the linear-phase FIR filter is obtained 
from the resultant approximated amplitude characteristic. 

Fig. 6 is a flowchart of the Remez Exchange algorithms 
while taking into account a frequency response of the 
pre-filter according to the present invention. 
20 The Remez Exchange algorithms while taking into 

account the concrete frequency response of the pre-filter 
is as follows : 
Step 0 

As illustrated in Fig. 6, firstly, an initial setting 



50 

is executed (F 101) . In this initial setting, setting of 
the linear-phase FIR filter, setting of a band, setting 
of a coefficient of the pre-filter, and setting of an 
initial extreme point are carried out. 

Items set concretely are as follows: 

• the number of tap , 

• even symmetry or odd symmetry of the linear-phase 
FIR filter, 

• the number of the band, 

• the frequency of both ends of the respective bands , 

• the desired amplitude value of the respective 
bands , 

• weighting to the respective bands, 

- coefficient of the pre-filter, 

- the frequency co (0) = co Jt t0) (k = 0 , .. , R) that 
becomes the extreme value in the approximation band. 

Note, the letter of right shoulder (i) expresses 
the number of repetition. 
Step 1 

Next, LaGrange' s interpolation polynomial equation is 
generated for interpolating an amplitude characteristic 
from the current extreme value point (F 102) . 

Above equation (5) indicates the target function of 
the Chebyshev approximation. The necessary and sufficient 
condition for obtaining the minimum value of the target 
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function of the Chebyshev approximation is indicated by 
the alternating theorem. Accordingly, parameter p (n) of 
the following equation is found in such a way that, on 
the basis of the alternating theorem, the weighted 
5 approximation error 6 (i) from the desired amplitude 

characteristic at the respective frequency points are of 
equal , and sign is alternated . 



j =3 p{e j< ")= ^p{n)co{tiCD) (12) 

Namely, the weighted approximation error of the 
□ 10 equation (7) at the frequency points co {i) = co k u) (k = 0, 

M: , R) satisfies the following equation. 

1 4"'f'^")-4 rf I-(- i )' 5< ' ) <*-°a-*) w 

15 Hereinafter, the right shoulder letter (i) is omitted 

for simplification. The equation (13) is carried as 
follows . 



20 



| P (,,M«oO + |^-^ M ) (k-0X.-,R) (14) 
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Matrix expression of the equation (14) is as follows 



1 cos(<w 0 ) cos(2cu 0 ) 
1 cos((0,) cos(2cy,) 

1 cos(w fl .,) cos(2cu /7 _,) 



cos((/? - l)a> 0 ) 
cos((/?- 1)^,) 



1 



W 



cosft/? -!)^.,) ^ jL, j 



(- 

1 cos(o; fl ) 005(2^) ••• cos((/?-lK) j 



p(o) 

p(/? " 1) 

<5 



(15) 



3(e M ) 

d(e J,,J ''-' ) 
D{e J '" R ) 



However, there are lot of amount of calculation for 
calculating this equation (15) , accordingly, 5 is 
analytically found firstly. 



8 = 



|(-iya,/#( e ^) 



(16) 



10 a *-rfe^n 



(17) 



(18) 
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This a, is surplus factor of component of k-column (R 

+ l,-row of the matrix F. Note, and t*. 

equation (8) and the equation (9) . 

Next, the following equation is provided by using 5 . 

C4 _^M)_(_ iy ^ (/C=0,...,*-1) (19) 

in order to find the amplitude characteristic of the 
frequency except for the extreme value point, the 
LaGrange's interpolation polynomial expression is used a 
10 an interpolation polynomial expression for interpolating 
using the extreme value point. Namely, A?' ) xs 
calculated in such a way that the interpolation is 
executed so as to obtain value C k at o» k (k=0 , . . , R-D 
by using the LaGrange' s interpolation polynomial 
15 expression. 



(20) 



x-x t 



x = coipj) 



(21) 
(22) 
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This result means the answer of the equation (15) , 
and satisfies the equation (13) automatically at the 
point of cor. 

Step 2 

5 The new extreme value point is determined from the 

amplitude characteristic obtained from the interpolation 
polynomial expression (F 103) , before returning to 
processing of F 102. Judgement is executed repeatedly 
whether or not the optimum approximation is obtained 
10 after position of the extreme value is not changed (F 
104) . 

Above described each extreme value point cak as a 
result of Step 1 is not necessarily the extreme value of 
the weighted function E{e JCO ) , in some cases, the point 
15 which becomes \E{e Jw )| > <5 (/) exists. Thus, the new extreme 
value a/' + l ' is determined from the all point simultaneous 
replacement method. 

All Point Simultaneous Replacement Method 
On the basis of the following equation, the extreme 
20 value of weighted approximation error calculated from the 
extreme value used for the interpolation is searched over 
the whole approximation bands , the searched extreme 
values are taken to be the new extreme values 

(,- + i) _ (/+i) (k = 0, 1, .. , R) f thus returning to Step 

CO - (O k 
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(23) 



5 It is decided that when the position of the extreme 

value has not been changed, the optimum approximation is 
obtained. This is the condition of the termination of the 
repetition, when the condition is satisfied, processing 
advances to next Step 3. 
10 Figs . 7A to 7C are conceptual views of the all point 

replacement method. 

Explanation is made simply. Black dots in Figs . 7A to 
7C explesse the extreme value points used for the 
interpolation. The weighted approximation error E^ 1 ) 

15 obtained from the extreme value points corresponds to the 
solid line. 

As illustrated in Fig.7A, values of the weighted 
approximation error at the extreme value points of the 
black dots become open circles. However, the actual 
20 extreme values are frequency indicated by black squares . 
Accordingly, the frequency indicated by the black squares 
is taken to be the new extreme value point, then 
returning to processing of Step 1 . 

Further, as illustrated in Fig.7B, since there is a 
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deviation of the frequency exists between the extreme 
value points used for the interpolation and the actual 
extreme value, the frequency indicated by the square is 
taken to be the new extreme value points , then returning 
to processing of Step 1. 

Then, as illustrated in Fig.7C, when the extreme 
value points used for the interpolation and the extreme 
value point (open circles) of the actual weighted 
approximation error become the same, and the repetition 
processing terminals . 
Step 3 

Coefficients of the linear-phase FIR filter are found 
from the approximated amplitude characteristic (F 105) . 
When n-tap of impulse response h (n) is found from 



the optimum approximation function r\e ), the n-tap of 
impulse response is obtained from the following equation 
instead of P (n) . 





Case 1 : 




Case 2 : 
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(26) 



Case 3 



/V-l 



*>--*S<-*^M^H)) 



(27) 



Case 4 



Further, when the impulse response k (n) of the 
10 equalizer of L = N - (U - 1) tap is found, the following 
equation is calculated. 



(29) 



15 



Case 1 : 



k(n) = K 
v ' L 



L-l 



(30) 



Case 2 : 
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(31) 



Case 3 : 

L-\ 




(33) 



Q As to the pre-filter, when the transfer function Z 

□ 10 (n) is 1 as indicated in the following equation, the 

H Remez Exchange algorithms in question is substantially 

equal to a normal Remez Exchange algorithms . 

Z(z) = l (34) 

15 

Here, comparison between the conventional technique 
and the method of the present invention is carried out by 
taking zero point for avoiding the Chessboard distortion 
to be the pre-filter . 
20 The frequency response of the pre-filter is expressed 

as follows : 



59 



IT! 



z( e jl " ) = (l + e- J '" +e- 2im +■■■ + e - Mu - l) ) 



1+ \ 2cos(mcL>), (7 ^ 35 ^ 

V 2cos|m + U : even 



The specification of the conventional technique and 
5 the present invention are indicated below. 

Fig. 8 is a view illustrating frequency response of 
low-pass filter which avoids the chessboard distortion 
designed depending on the present invention. Figs . 9A and 
9B are views for comparing the frequency response of the 
10 low-pass filter designed depending on the conventional 

technique with the present invention respectively. Fig. 10 
is a view illustrating weighted approximation error to 
the filter designed depending on the present invention. 
Specification of the conventional technique 
15 Linear-phase FIR filter of equalizer 

• 22 taps 

• even symmetry 
Method of Design 

The equalizer is made to design with the Remez 
20 Exchange algorithms, thus the pre-filter of zero point is 
folded into the linear-phase FIR filter in order to avoid 
the chessboard distortion of U = 3. Finally, the linear- 
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phase FIR filter with 24 taps is designed. 



[Table 1] 

Band 



Band 


Range of frequency 


gain 


weight 


1 


0<co^0 . 3n 


1 


1 


2 


O.Sn^co^n 


0 


1 



In Figs . 9A and 9B, the conventional technique is 
illustrated by dotted line written as "Conventional". 
Specification of the present invention 

Linear-phase FIR filter 

- 2 4 taps 

• even symmetry 

■ U = 3 

Method of Design 

Designing has been attempted with Remez Exchange 
algorithms by taking into account the frequency respon 
of the pre-filter of zero point for avoiding the 
chessboard distortion. 



61 

[Table 2] 



Band 



Band 


Range of frequency 


gain 


weight 


1 


0^(0^0 . 3n 


1 


1 


2 


0 . Sn^co^n 


0 


1 



The dotted lines in Fig. 8 and Figs. 9A and 9B 
indicate the frequency characteristic of the pre-f ilter 
10 and the frequency characteristics of the equalizer of the 
zero point for avoiding the chessboard distortion. The 
solid lines indicate both the frequency characteristic 
(Proposed H (z)) obtained finally and the frequency (Zero 
Point) that must be H (z) = 0 for avoiding the chessboard 
15 distortion. 

As seen from Fig . 8 to Fig. 10, the method of the 
present invention is capable of obtaining good frequency 
response characteristic compared with the conventional 
technique . 

20 Next, as a modified example, an algorithm which 

satisfies an attenuation quantity of the designated stop 
band will be explained. 

This algorithms is an algorithm satisfies the 
attenuation quantity (weighted approximation error) of 

25 the designated stop band in connection with the low-pass 
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filter with 2 numbers of bands or the high-pass filter by 
using the Remez Exchange algorithms which takes into 
account the frequency response of the pre-filter. 

Further, in later explanation, description is made 
5 for the low-pass filter, however, when the algorithms is 
intended to apply to the high -pass filter, the pass -band 
and the block-band are inversely related. 

As a method for satisfying an attenuation quantity of 
the stop band, there exists three kinds of approaches. 
10 The first approach is that the end frequency cop of 

the pass band is made variable, the starting frequency cos 
of the stop band is made to be fixed, and the number of 
taps is made to be fixed. 

The second approach is that the end frequency cop of 
15 the pass band is made to.be fixed, the starting frequency 
cos of the stop band is made variable, and the number of 
taps is made to be fixed. 

The third approach is that the end frequency cop of 
the pass band is made to be fixed, the starting frequency 
20 cos of the stop band is made to be fixed, and the number 
of taps is made variable. 

The first, the second, and the third approaches are 
explained below with reference to the accompanying 
drawings successively. 
25 First Approach 
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In the first approach, the maximum end frequency cop 
of the pass band that satisfies the attenuation quantity 
dBs of the designated stop band is found. 

Fig. 11 is a view illustrating flowchart of algorithms 
5 for finding filter that satisfies the attenuation 
quantity of the stop band. 

Fig. 12 is a view illustrating flowchart of algorithms 
for finding filter that has the maximum end frequency of 
the pass band, which end frequency satisfies the 
10 attenuation quantity of the stop band, in which free 

parameters (variable) and fixed parameter (variable) are 
expressed . 

Here, free parameters in the algorithms, object, and 
the principle of the algorithms are enumerated as 
15 follows : 

* Free parameters: end frequency cop of the pass band. 

* Object: it is an object to obtain a filter that has 
the maximum end frequency cop of the pass band, which end 
frequency satisfies the attenuation quantity of the 

20 designated stop band. 

* Principle: the starting frequency of the pass band 
and the frequencies of both ends of the stop band are 
fixed, and the end frequency of the pass band is the free 
parameter . In the Chebyshev approximation according to 

25 the Remez Exchange algorithms , 
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• the end frequency cop of the pass band moves far 
away rather than the starting frequency us 

-►the attenuation quantity of the stop band becomes 

large . 

5 • the end frequency cop of the pass band approaches 

the starting frequency cos of the stop band 

the attenuation quantity of the stop band 
becomes small . 

Namely, co p (cur) is made to prepare at the frequency (in 
10 the vicinit-y of co = 0) which is far from the starting 
frequency cos of the stop band. While, co p <pre> which is 
close to the starting frequency of the stop band is made 
to prepare as initial frequency. Then, the position cop of 
the maximum end frequency is found within the pass band 
15 that satisfies the designated attenuation quantity while 
utilizing the dichotomizing method. 

Further, the golden section method is the most 
efficient method in the line search methods of such 
parameter. However, the dichotomizing (two branching) 
20 method is adopted here. The dichotomizing method has the 
algorithms which is easy to understand. 

Fig. 11 and contents of each step processing F 102, F 
103, F 104, F 105 of the algorithms described hereinafter 
are the same as the Remez Exchange algorithms under 
25 taking into account the frequency response of the pre- 
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filter described in relation to Fig . 6 . Accordingly, the 
same signs as that of Fig. 6 are used in connection with 
these processing. 
Step 10 

As illustrated in Fig. 11, firstly, initial setting is 
executed (F 201) . In this initial setting, setting of the 
linear-phase FIR filter, setting of the band, setting of 
the coefficient of the pre-filter, setting of the initial 
extreme value point, designation of the attenuation 
quantity of the stop band, and setting of the initial 
frequency of the dichotomizing method are executed. 

Concrete items thus set are as follows: 

• the number of tap 

• linear-phase FIR filter with even symmetry or odd 
symmetry 

• the number of the band is two 

• the starting frequency of the pass band 

• gain of the pass band 

• frequencies of both ends of the stop band 

• gain of the stop band 

• weighting in connection with the pass band and the 
stop band 

• coefficients of the pre-filter 

• attenuation quantity dBs (namely, pointing 
largeness 52 of ripple of the stop band) 
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• frequency u <0> = a> k t0) (k = 0 , . . , R) that becomes 
the extreme value in the approximation band. Provided 
letter of right shoulder (i) expresses the number of 
repetition. 

Also, Fig. 13 is a view illustrating initial frequency 
of the dichotomizing method in the algorithms for finding 
filter having the maximum end frequency of the pass band 
that satisfies the attenuation quantity of the stop band. 

As illustrated in Fig. 13, in the present example, the 
following values are given as the initial frequencies of 
the dichotomizing method. 

O^'^[o] = the v i c i n ity of starting frequency of the 
stop band (37) 

a^"")[o] = the vicinity of starting frequency of the 
pass band (38) 

Portion of [t] at the back of up expresses the number 
of cycle. 

Here, the Remez Exchange algorithms taking into 
account the frequency response of the pre-filter is 



the attenuation quantity dB pre [OJ in the stop band is 
found. The following steps are explained on the basis of 



executed to the frequency 




thus the minimum of 
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the assumption. 
Ster> 11 

The Remez Exchange algorithms taking into account the 
frequency response of the pre-filter is executed (F 102, 
5 F 103, F 104) . 

Concretely, in the processing F 102, the 
interpolation polynomial expression is generated for 
interpolating an the amplitude characteristic from the 
extreme value point at the time of a>^ £/r [f] . 
10 Subsequently , in the processing F 103, the new 

extreme value point is determined from the amplitude 
characteristic obtained from the interpolation polynomial 
expression . 

Then, in the processing F 104, the repetition 
15 judgement of the Remez Exchange algorithms is executed. 
Step 12 

Next, the attenuation quantity of the stop band is 
examined (F 20 6) . 

The minimum attenuation quantity (the maximum 
20 weighted approximation error 52) <^£^[f] is examined in 
the stop band by using the interpolation polynomial 
expression founded in the processing F 102. 
Step 13 

Next, comparison is attempted in connection with the 
25 attenuation quantity of the designated stop band (F 207) . 
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Concretely, in the case that the result of the 
comparison satisfies the following equations (39) or (40) 
while comparing with the attenuation quantity dBs of the 
designated stop band, processing is moved to processing 
5 of Step 15 (F 105) . While when the result of the 

comparison does not satisfy equations (39) or (40) , 
processing is moved to processing of Step 14 (F 208) . 

dBf u) \t\-dB s \{e x (39) 
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Up 



i^'W-OI^Wl^ (40) 



Provided that el and e2 are extremely small values. 
Step 14 

Setting of the band is changed when the result of the 
15 comparison does not satisfy equations (39) or (40) , in 
the comparison to the attenuation quantity dBs of the 
designated stop band (F 208) . 

Concretely, the end frequency ^^[t + l] of new pass 
band is made to set. Explanation is made while dividing 
20 into two cases of the first loop and the second loop as 
the setting method. 
First time: 

In the case of the first time, there are three cases 
illustrated in Figs.l4A to 14C. 



IS 
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Namely, cases illustrated in Fig.l4A, Fig.l4B, and 
Fig.l4C are considered in connection with the frequencies 

The case illustrated in Fig.l4A is one in which the 
5 both satisfy the designated attenuation quantity. In th 

case, the frequency G^/"^[o] close to the starting 
frequency cos of the stop band is taken to be solution, 
thus advancing toward to the processing of Step 15. 

The case illustrated in Fig.l4B is one in which the 
10 both do not satisfy the designated attenuation quantity. 
In this case, the present number of taps can not realize 
the designated attenuation quantity. Accordingly, that 
effect is indicated to end processing. 

The casze illustrated in Fig.l4C is one in which only 
15 one side satisfies the designated attenuation quantity. 
In this case, processing is moved to Step 11 in such a 
way as following equation (41) . 

Further, in the present method for giving initial 

frequency point, there is no case in which Ctfjf ^[o] is 

20 satisfied and <^T''[P] is not satisfied. 

cujr = 0.5 x (cu / ( r ) [0]+ ^[oj (41) 
dB lgw) [l]^dB [cur) [o] 
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Second time and later : 



In the case of second time and later, two cases 
indicated in Figs . 15A and 15B are considered. In the 
determining method of the new frequency in the second 
5 time and later, the frequency which satisfies the 
designated attenuation quantity dBs is preserved 

invariably within CO^ ie ^\t +l] . The case illustrated in 
Fig.lSA, Fig.lSB is considered in connection with the 



side satisfies the designated attenuation quantity. In 
this case, the processing is moved to next step 11 in 
such a way as the following equation (42) . 

Further, since the frequency which satisfies the 
15 designated attenuation quantity always conforms to 



frequencies CD 



t~>[ t ] w <~)[ t ]. 
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The case illustrated in Fig.lSA is one in which one 




there does 0) 




not exist the case that does 



not satisfy. 




(42) 



dB {prc) [t+l] = dB pre [t] 



20 



The case illustrate in Fig.lSB is one in which the 



both satisfies the designated attenuation quantity. In 
this case, the processing is moved to next step 11 in 
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such a way as the following equation (43) 



dB {pre) [t + 1] = dB [cur) [t] 



5 Step 15 

The coefficients of the linear-phase FIR filter is 
obtained from the approximated amplitude characteristic 
(F 105) . 

Namely, the coefficients of the linear-phase FIR 
10 filter is obtained from the amplitude characteristic 
obtained finally. 

Fig. 16 is a view illustrating frequency response 
characteristic of the low-pass filter obtained from the 
algorithms for obtaining "filter having the maximum end 
15 frequency of the pass band that satisfies the attenuation 
quantity of the stop band" while taking the filter of the 
zero point for avoiding the chessboard distortion. 

Fundamental algorithms , variable intended to find, 
and band of this case are as follows : 
20 Fundamental algorithms 

It is the Remez Exchange algorithms taking into 
account frequency characteristic of the pre-f ilter . 
• 2 4 taps 
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• even symmetry 

• U = 3 (filter coefficient is executed U times 
such that DC gain becomes U) 

• attenuation quantity of the stop band is not more 
5 than - 40 dB 

Variable intended to find 
It is the end frequency cop of the pass band. 

[Table 3] 



Band 



band 


Frequency region 


gain 


weight 


Pass band 


0<<o£co p 


i 


1 


Stop band 


0 . Sn^co^n 


0 


1 



Further, in Fig. 16, the solid line indicates 
15 frequency response of low-pass filter that has the 

maximum end frequency of the pass band which satisfies 
the attenuation quantity of the stop band. Furthermore, 
the dotted line indicated by (D indicates frequency which 
becomes zero point for avoiding the chessboard 
20 distortion, and the dotted line indicated by <D 
indicates period of the band given beforehand. 

As seen from Fig. 16, the low-pass filter which is 
obtained from the algorithms for finding "filter having 
the maximum end frequency of the pass band which 
25 satisfies the attenuation quantity of the stop band" 
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according to the present invention has a good frequency 
response characteristic . 
Second approach 
Next, in the second approach, the minimum starting 
5 frequency cos of the stop band that satisfies the 

attenuation quantity dBs of the designated stop band is 
found. 

The flowchart for this attempt is equivalent to that 
of Fig. 11 referred to in the first approach. 
10 Fig. 17 shows a free parameter (variable) and fixed 

parameter (variable) in the algorithms for finding filter 
having the minimum starting frequency of the stop band 
that satisfies the attenuation quantity of the stop band. 
Here, the free parameter of this algorithms , objects, 
15 and the principle of the. algorithms thus enumerated are 
as follows : 

* free parameter: it is the starting frequency cos of 
the stop band. 

* objects: it is an object to obtain filter having 
20 the minimum starting frequency cos of the stop band that 

satisfies attenuation quantity of the designated stop 
band . 

* principle: the end frequency of the stop band and 
the frequencies of both ends of the pass band are fixed, 

25 and the starting frequency of the stop band is the free 
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parameter . In the Chebyshev approximation according to 
the Remez Exchange algorithms, 

• the starting frequency cos of the stop band moves 
far away rather than the end frequency cop of the pass 

5 band -> the attenuation quantity of the stop band becomes 
large . 

• the starting frequency cos of the stop band 
approaches the end frequency cop of the pass band — ► the 
attenuation quantity of the stop band becomes small. 

10 Namely, co g (pi:e) is made to prepare that is far from the 

end frequency of the pass band. While, co a (cur) which is 
close to the end frequency of the pass band is made to 
prepare as initial frequencies. Then, position (as of the 
minimum starting frequency is found within the stop band 

15 that satisfies the designated attenuation quantity while 
utilizing the dichotomizing method. 

Further, the golden section method is the most 
efficient method in the line search methods of such 
parameter. However, the dichotomizing method is made to 

20 adopt here. The dichotomizing method has the algorithms 
which is easy to understand. 

Also, Fig. 11 and contents of each step processing F 
102, F 103, F 104, F 105 of the algorithms described 
hereinafter are the same as the Remez Exchange algorithms 

25 while taking into account the frequency response of the 
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pre-filter described in relation to Fig. 6 in the same way 
as the case of the first approach. Accordingly, the same 
signs as that of Fig . 6 are used in connection with these 
processing . 
5 Step 20 

As illustrated in Fig. 11, firstly, initial setting is 
executed (F 201) . In this initial setting, setting of the 
linear-phase FIR filter, setting of the band, setting of 
^ a coefficient of J the pre-filter, setting of the initial 

; -_ s 

g 10 extreme value point, designation of attenuation quantity 

;n of the stop band, and setting of initial frequency of 

dichotomizing method are executed. 
;m Concrete items thus set are as follows: 

£J * the number of taps 

:=k * 15 • linear-phase FIR filter with even symmetry or odd 

symmetry 

• the number of the bands is two 

• the starting frequency of the pass band 

• gain of the pass band 

20 • frequencies of both ends of the stop band 

• gain of the stop band 

• weighting to the pass band and the stop band 

• coefficients of the pre-filter 

• attenuation quantity dBs (namely, pointing 
25 largeness 52 of ripple of the stop band) 
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• frequency G) (0) = ca k <0) (k = 0, . . , R) that becomes 
the extreme value in the approximation band. Provided 
letter of right shoulder (i) expresses the number of 
repetition . 

5 Also, Fig. 18 is a view illustrating initial frequency 

of the dichotomizing method in the algorithms for finding 
filter having the minimum starting frequency of the stop 
band that satisfies the attenuation quantity of the stop 
band . 

10 As illustrated in Fig. 18 , in the present example, the 

following values are given as the initial frequencies of 
the dichotomizing method. 

the vicinity of end frequency of the pass 
15 band (44) 

co[ cw/ '[o]= the vicinity of end frequency of the 
stop band (45) 

Portion of [t] at the back of cos expresses the number 
20 of cycle . 

Here, the Remez Exchange algorithms taking into 
account the frequency response of the pre-filter is 
executed to the frequency Cl> j pn? ^[o] , thus the minimum of 
the attenuation quantity a^ pn2 ^[o] in the stop band is 
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found. The following steps are explained on the basis of 
this assumption . 
Step 21 

The Remez Exchange algorithms taking into account the 
5 frequency response of the pre-filter is executed (F 102, 
F 103, F 104) . 

Concretely, in the processing F 102, the 
interpolation polynomial expression is generated for 
interpolating an amplitude characteristic from the 



10 extreme value point at the time of (D s \y\ . 

Subsequently , in the processing F 103, the new 
extreme value point is determined from the amplitude 
characteristic obtained from the interpolation polynomial 
expression. 

15 Then, in the processing F 104, repetition judgement 

of the Remez Exchange algorithms is executed. 
Step 22 

Next, the attenuation quantity of the stop band is 
examined (F 206) . 
20 The minimum attenuation quantity (the maximum 



the stop band while using the interpolation polynomial 



GO 




weighted approximation error 52) 




expression founded in the processing F 102. 



Step 23 



78 

Next, comparison to the attenuation quantity of the 
designated stop band is executed (F 207) . 

Concretely, in the case that result of the comparison 
satisfies following equations (46) or (47) while 
5 comparing with the attenuation quantity dBs of the 

designated stop band, processing is moved to processing 
of Step 25 (F 105) . While when result of the comparison 
does not satisfy equations (46) or (47) , processing is 
moved to processing of Step 24 (F 208) . 
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d$ a, *\t]-dB s \(e x (46) 



o>]r\ty to^\t%B t . (47) 



15 



Provided that el and e2 are extremely small values. 
Step 24 

Setting of the band is changed when the result of the 
comparison does not satisfy equations (46) or (47) , in 
the comparison to the attenuation quantity dBs of the 
20 designated stop band (F 208) . 

Concretely, the starting frequency a^ CMr *[f+l] of new 
stop band is made to set. Explanation is made while 
dividing into two cases of the first loop and the second 
loop as the setting method- 
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First time: 

In the case of the first time, there are three cases 
illustrated in Figs.l9A to 19C. 

Namely, cases illustrated in Fig.l9A, Fig.l9B, and 



The case illustrated in Fig.l9A is one in which the 
both satisfy the designated attenuation quantity. In this 
case, the frequency 0>j pn ^[o] close to co=0 is taken to be 
solution, thus advancing toward to the processing of Step 
25 . 

The case illustrated in Fig.l9B is one in which the 
both do not satisfy the designated attenuation quantity. 
In this case, the present number of taps can not realize 
the designated attenuation quantity. Accordingly, that 
effect is indicated to end. 

The case illustrated in Fig.l9C is one in which only 
one side satisfies the designated attenuation quantity. 
In this case, processing is moved to Step 21 in such a 
way as following equation (48) . 

Further, in the present method for giving initial 



Fig.l9C are considered to frequencies 




frequency point, the case in which 




satisfied 



and CO 




is not satisfied does not exist. 
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aj { r' ] [\. | = 0.5 x (of ttr, [o]+ o^o) (48) 
dBl /frc) [l] = dBr[0] 

Second time and later : 
In the case of second time and later, two cases 
5 indicated in Figs.20A and 20B are considered. In the 
determining method of the new frequency in the second 
time and later , the frequency which satisfies the 
attenuation quantity dBs thus designated is preserved 

invariably within 0^ p/ ^[f+l]. The case illustrated in 
10 Fig.20A, Fig.20B is considered to the frequencies 




The case illustrated in Fig.20A is one in which one 
side satisfies the designated attenuation quantity. In 
this case, the processing is moved to next step 21 in 
15 such a way as the following equation (49) . 

Further, since- the frequency which satisfies the 
designated attenuation quantity always conforms to 

GJ^[*], the case that W^l/] does not satisfy does not 
exist . 



col'" c] [i + l]=(ol" re) [t] 

aj [ r ] [( + l] = 0.5 x (a>i c " r) [f ]+ (o ( r c) [t^ < 49 ) 

clB ( J" a) [[ +l] = dB , f" r) [t] 
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The case illustrate in Fig. 2 OB is one in which the 
both satisfies the designated attenuation quantity. In 
this case, the processing is moved to next step 21 in 
such a way as the following equation (50) . 

oj^t + 1] = co ( s cur) [t] 

0) { s CUr) [t + 1] = <D { s CUr) [t] + 0.5 X [w { 5 CUr) [t\- (D { s Pr * ] [t]) (50) 

dB [ s pre) [t + 1] = dB [ 5 cur \t] 

Step 25 

The coefficients of the linear-phase FIR filter are 
found from the approximated amplitude characteristic (F 
105) . 

Namely, the coefficients of the linear-phase FIR 
filter are found from the amplitude characteristic 
obtained finally. 

Fig. 21 is a view illustrating frequency response 
characteristic of the low-pass filter obtained from the 
algorithms for finding "filter having the minimum 
starting frequency of stop band which filter satisfies 
the attenuation quantity of the stop band" while taking 
the filter of the zero point for avoiding the chessboard 
distortion . 

Fundamental algorithms, variables intended to find, 
and band of this case are as follows : 



10 
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Fundamental algorithms 
It is the Remez Exchange algorithms taking into 
account frequency characteristic of the pre-f ilter . 

• 2 4 taps 

• even symmetry 

• U = 3 (filter coefficient is executed U times 
such that DC gain becomes U) 

• attenuation quantity of the stop band is not more 
than - 40 dB 

Variable intended to find 
It is the starting frequency ws of the stop band. 



15 



[Table 4] 



Band 



band 


Frequency region 


gain 


weight 


Pass band 


0<co^0 . 3n 


1 


1 


Stop band 




0 


1 



Further, in Fig. 21, the solid line indicates 
20 frequency response of low-pass filter that has the 
minimum starting frequency of the stop band which 
satisfies the attenuation quantity of the stop band. 
Furthermore, the dotted line indicated by (D indicates 
frequency which becomes zero point for avoiding the 
25 chessboard distortion, and the dotted line indicated by 



83 

® Indicates period of the band given beforehand. 

As seen from Fig. 21, the low-pass filter which is 
obtained from the algorithms for finding "filter having 
the minimum starting frequency of the stop band which 
5 satisfies attenuation quantity of the stop band" 

according to the present invention has a good frequency 
response characteristic. 
Third approach 
Finally, in the third approach, the minimum number N 
10 of taps that satisfies the attenuation quantity dBs of 
the designated stop band is found. 

Fig. 22 is a view illustrating flowchart of the 
algorithms for finding filter with the minimum number of 
tap that satisfies the attenuation quantity of the stop 
15 band. 

Here, free parameters of this algorithms, objects, 
and the principle of the algorithms thus enumerated are 
as follows : 

* free parameters: it is the number of taps 

20 * objects: it is an object to obtain filter having 

the minimum number of tap that satisfies attenuation 
quantity of the designated stop band. 

* principle: since variable of the band is fixed, 
when it is impossible to satisfy the attenuation quantity 

25 of the designated stop band, the number of tap is 
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increased by one tap . 

Also, Fig. 22 and contents of each step processing F 
102, F 103, F 104, F 105 of the algorithms described 
hereinafter are the same as the Remez Exchange algorithms 
5 while taking into account the frequency response of the 
pre-filter described in relation to Fig. 6 in the same way 
as the case of the first approach. Further, contents of 
processing F 206 is the same as the processing explained 
in relation to Fig. 11 in the same way as the case of the 
10 second approach. Accordingly, the same signs as that of 
Fig. 6 and Fig . 11 are used in connection with these 
processing . 

Step 30 

As illustrated in Fig. 22, firstly, initial setting is 
15 executed (F 401) . In this initial setting, setting of 
linear-phase FIR filter, setting of band, setting of 
coefficient of the pre-filter, setting of the initial 
extreme value point, and designation of attenuation 
quantity of the stop band, are executed. 
20 Concrete items thus set are as follows : 

• the number of initial tap 

• linear-phase FIR filter with even symmetry or odd 
symmetry 

• the number of the band is two 

25 • the starting frequency and end frequency of 



respective bands 

• gain of respective bands 

• weighting of respective bands 

• coefficient of the pre-filter 

5 • attenuation quantity dBs (namely, pointing 

largeness 52 of ripple of the stop band) 

• frequency g> (0> = w k <0) (k = 0 , . . , R) that becomes 
Jj the extreme value in the approximation band. Provided 
y letter of right shoulder (i) expresses the number of 

Q 10 repetition. 

^ Step 31 

J J The Remez Exchange algorithms taking into account the 

frequency response of the pre-filter is executed (F 102, 

J F 103, F 104) . 

15 Concretely, in the processing F 102, the 

interpolation polynomial expression is generated for 
interpolating an amplitude characteristic from the 

(cur)r.l 

extreme value point at the time of C0 p [tj . 

Subsequently, in the processing F 103, the new 
20 extreme value point is determined from the amplitude 

characteristic obtained from the interpolation polynomial 
expression . 

Then, in the processing F 104, repetition judgement 
of the Remez Exchange algorithms is executed. 
25 Step 32 
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Next, the attenuation quantity of the stop band is 
examined (F 206) . 

The minimum attenuation quantity (the maximum 
weighted approximation error 52) </£?j clir) [f] is examined in 
5 the stop band while using the interpolation polynomial 
expression founded in the processing F 102. 
Step 33 

Next, comparison to the attenuation quantity of the 
designated stop band is executed (F 412) . 

10 Concretely, in the case that result of the comparison 

satisfies the following equations (51) while comparing 
with the attenuation quantity dBs of the designated stop 
band, processing is moved to processing of Step 35 (F 
105) . While when result of the comparison does not 

15 satisfy equations (51) , processing is moved to processing 
of Step 34 (F 413) . 

\d4 cu ^dB s \ (51) 

Step 34 

20 One tap is increased (F 413) . 

Namely, one tap is increased in addition to the 
present number of tap, thus processing is moved to 
processing of Step 31, 
Step 35 

25 The coefficients of the linear-phase FIR filter is 
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found from the approximated amplitude characteristic (F 
105) . 

Namely, the coefficients of the linear-phase FIR 
filter is found from the amplitude characteristic 
5 obtained finally. 

Fig. 23 is a view illustrating frequency response 
characteristic of the low-pass filter obtained from the 
algorithms for finding "the minimum number of tap that 
satisfies attenuation quantity of the stop band" while 
10 taking the filter of the zero point for avoiding the 
chessboard distortion. 

Fundamental algorithms, variables intended to find, 
and band of this case are as follows : 
Fundamental algorithms 
15 it is the Remez Exchange algorithms taking into 

account frequency characteristic of the pre-f ilter . 

• the number of initial tap is 10 taps 

• even symmetry 

• U = 3 (filter coefficient is executed U times 
20 such that DC gain becomes U) 

• attenuation quantity of the stop band is not more 
than - 60 dB 

Variables intended to find 
It is the number N of tap. 

25 
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[Table 5] 

Band 



band 


Frequency region 


gain 


weight 


Pass band 


O^co^O . 3n 


1 


1 


Stop band 


0 . Sn^co^n 


0 


1 



Further, in Fig. 23, the solid line indicates 
frequency response of low-pass filter of the minimum 
; ! 2 number of tap (36 taps) in which the attenuation quantity 

10 of the stop band becomes not more than -60 dB . 
S Furthermore, the dotted line indicated by <D indicates 

frequency which becomes zero point for avoiding the 
□ chessboard distortion, and the dotted line indicated by 

Q (g) indicates period of the band given beforehand, 

p 15 As known from Fig. 23, the low-pass filter which is 

obtained from the algorithms for finding "the minimum 
number of tap which realizes attenuation quantity of the 
stop band" according to the present invention has good 
frequency response characteristic. 
20 Subsequently, as the second modified example, 

explanation is performed about algorithms for finding 
filter that causes the signal to pass frequency point of 
transition band that satisfies attenuation quantity of 
the stop band. 
25 Here, followings are free parameter in the 
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algorithms, objet, and the principle of the algorithms. 

* Free parameter: end frequency cap of the pass band 
and starting frequency cos of the stop band. 

* Object: it is an object to determine a band that 

5 satisfies attenuation quantity dBs and makes attenuation 
quantity dBc with frequency coc of transition band pass 
through. Namely, it is an object to obtain the maximum 
end frequency cop of the pass band and the minimum 
starting frequency cos of the stop band by which the 
10 attenuation quantity dBc is realized with specific 
frequency coc of the transition band. 

* Principle: the starting frequency of the pass band 
and the end frequency of the stop band are fixed, and the 
end frequency cop of the pass band and the starting 

15 frequency cos of the stop band are the free parameters. 

Since there are two free parameters, it is impossible to 
determine appropriately if moving simultaneously. 
Accordingly, one side of parameter is fixed, while the 
other side of parameter that satisfies attenuation 
20 quantity of the stop band is found. When the band does 

not make the attenuation quantity dBc with the frequency 
coc of the transition band pass through, it causes the 
fixed parameter to be changed. Above repetition is 
carried out so that the band which makes the attenuation 
25 quantity dBc with the frequency coc of the transition band 
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pass through is determined. 

In the method for finding parameter that satisfies 
the attenuation quantity of the stop band, there are two 
kinds of approaches . 
5 Namely, the first approach is that filter which has 

the minimum starting frequency of the stop band that 
satisfies the attenuation quantity of the stop band is 
found . 

The second approach is that filter which has the 
10 maximum end frequency of the pass band that satisfies the 
attenuation quantity of the stop band is found. 

Algorithms on the basis of the above two kinds of 
approaches are explained successively. 

Fig. 2 4 is a flowchart of algorithms for finding 
15 filter 

That satisfies the attenuation quantity of the stop band 
and that enables the signal to be passed through the 
frequency point of the transition band. 

The algorithms on the basis of the filter having the 
20 minimum starting frequency of the stop band that 

satisfies the attenuation quantity of the stop band will 
be explained in relation to Fig. 24 to Fig. 29. 

Namely, the present algorithms utilizes the 
algorithms for finding "filter having the minimum 
25 starting frequency of the stop band that satisfies the 
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attenuation quantity of the stop band". 

The concrete course of the algorithms , as illustrated 
in Fig. 25, cos is obtained with aforementioned algorithms 
as inner loop, further loop for finding cop is put on 
5 outside thereof to determine parameter of the band. 

Namely, the outside loop causes the end frequency cop of 
the pass band to be fixed, while the inner loop causes 
the starting frequency cos of the stop band to be found 
that satisfies the attenuation quantity dBs of the stop 
10 band. 

When the amplitude characteristic thus found does not 
pass through the designated point (coc, dBc) of the 
transition band, cop is searched while using dichotomizing 
method like the case of searching "filter having the 

15 maximum end frequency of the pass band that satisfies the 
attenuation quantity of the stop band" . 

Hereinafter, concrete algorithms is indicated, 
however, since the algorithms for finding "filter having 
the minimum starting frequency of the stop band that 

20 satisfies the attenuation quantity of the stop band" to 
be the inner loop is the same as that described above, 
thus description thereof is omitted. 

Also, the content of Fig. 24 and each step processing 
F 102, F 103, F 104, and F 105 of the algorithms 

25 explained below, like the case of the first approach, is 
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the same as the Remez Exchange algorithms taking into 
account the frequency response of the pre-f ilter 
explained in relation to Fig. 6. Further, the content of 
each processing F 206, F 207, F 208, like the case of the 
5 second approach, is the same as processing explained in 
relation to Fig. 11 namely the algorithms of "filter 
having the minimum starting frequency of the stop band 
that satisfies the attenuation quantity of the stop 
band". Accordingly, as for this processing, the same 
10 marks as Fig. 6 and Fig. 11 are used. 
Step 40 

As illustrated in Fig. 24, firstly, initial setting is 
executed (F 301) . In this initial setting, setting of 
linear-phase FIR filter, setting of band, setting of 
15 coefficient of the pre-f ilter, setting of the initial 

extreme value point, designation of attenuation quantity 
of the stop band, and setting of initial frequency of 
dichotomizing method are executed. 

Concrete items thus set are as follows: 
20 • the number of tap 

• linear-phase FIR filter with even symmetry or odd 

symmetry 

• the number of the band is two 

• the starting frequency of the pass band « = 0 
25 • gain of the pass band 
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• the end frequency of the stop band to = n 



• gain of the stop band 



• weighting to the pass band and the stop band 



• coefficient of the pre-filter 



5 



attenuation quantity dBs of the stop band (namely, 



pointing largeness 52 of ripple of the stop band) 

* frequency coc of the transition band and attenuation 
quantity dJBc thereof 

• frequency co <0> = co k (0) (k = 0, , R) that becomes 
10 the extreme value in the approximation band. Provided 

letter of right shoulder (i) expresses the number of 
repetition. 



of the dichotomizing method in the algorithms for finding 
10 filter that causes the signal to pass the frequency point 

of the transition band, and that satisfies the 

attenuation quantity of the stop band. 

As illustrated in Fig. 26, in the present example, 

following values are given as the initial frequencies of 
20 the dichotomizing method. 



Also, Fig. 26 is a view illustrating initial frequency 




the vicinity of starting frequency of the 



pass band 



(52) 




the vicinity of coc 



(53) 
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Portion of [t] at the back of cop expresses the number 
of cycle. 

Here, the Remez Exchange algorithms taking into 
account the frequency response of the pre-filter is 
5 executed to the frequency &>p Pr ^[o] , thus the minimum 

starting frequency toj pr< ^[o] of the stop band, that 
satisfies the designated attenuation quantity dBs in the 
stop band is found, the gain fl49^ cur ^0] corresponding to coc 
at this time is found. The following steps are explained 
10 on the basis of this assumption. 
Step 41 

The algorithms for finding "filter having the minimum 
starting frequency of the stop band that satisfies the 
attenuation quantity of the stop band" is executed (F 

15 102, F 103, F 104, F206, F 207, F 208). 

Concretely, in the processing F 102, interpolation 
polynomial expression is generated for interpolating an 
amplitude characteristic from the extreme value point at 
the time of ^^[f] - 

20 Subsequently, in the processing F 103, the new 

extreme value point is determined from the amplitude 
characteristic obtained from the interpolation polynomial 
expression . 

Then, in the processing F 104, repetition judgement 
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of the Remez Exchange algorithms is executed. 

Next, in the processing F 206, the minimum 
attenuation quantity is found in the stop band. 

Next, in the processing F 207, end condition of the 
5 algorithms for finding filter having the minimum starting 
frequency cos of the stop band that satisfies the 
attenuation quantity of the stop band is obtained. 

Also, in the processing F 208, setting of the 
starting frequency Ctfj c< "^ of the new stop band is 
10 executed. 

Step 42 

Next, the attenuation quantity of the frequency that 
is designated as the transition band (F 309) is examined. 
In the processing F 102, the attenuation quantity 
15 fl^ cw/ ^[t] to the frequency wc which is designated as the 
transition band while using LaGrange' s interpolation 
polynomial expression found finally with the algorithms 
for finding "the minimum starting frequency of the stop 
band that satisfies the attenuation quantity of the stop 
20 band" is examined. 
Step 43 

Comparison to the designated attenuation quantity of 
the transition band is executed (F 310) . 

Concretely, in the case that result of the comparison 
25 satisfies following equations (54) or (55) processing is 
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moved to processing of Step 45 (F 105) . While when result 
of the comparison does not satisfy equations (54) or 
(55) , processing is moved to processing of Step 44 (F 
311) . 



\d4 tw) \f]-dB c \(e x (54) 



(55) 



Provided that el and e2 are extremely small values. 
Step 44 

Setting of the band is changed when result of the 
comparison does not satisfy equations (54) or (55) , in 
the comparison between the frequency coc of the transition 
band and the attenuation quantity dBc (F 311) . 

Concretely, the end frequency W^ ur \t + l] of new pass 
band is made to set. Explanation is performed while 
dividing into two cases of the first loop and the second 
loop as the setting method. 
First time: 

In the case of the first time, there are three cases 
illustrated in Figs.27A to 27C. 

Namely, cases illustrated in Fig.27A, Fig.27B, and 
Fig.27C are considered to frequencies 
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I" a 



15 



The case illustrated in Fig.27A is one in which the 
both satisfy the designated attenuation quantity dBc . In 



solution, thus advancing toward the processing of Step 
45 . 

The case illustrated in Fig.27B is one in which the 
both do not satisfy the designated attenuation quantity. 
In this case, the present number of taps can not realize 
the designated attenuation quantity. Accordingly, that 
effect is indicated to end. 

The case illustrated in Fig.27C is one in which only 
one side satisfies the designated attenuation quantity. 
In this case, processing is moved to Step 41 in such a 
way as following equation (56) . 

Further, in the present method for giving initial 
frequency point, there does not exist the case in which 
<^ Ct/r) [°] is satisfied and 0^ pre) [o] is not satisfied. 



this case, the large frecjuency (x) 



is taken to be 




w ;r ) [i]=o.5x( (y ;r' 0 [o] + a;;r ) [o]) 



(56) 



dB[ pre) [l\ = dB { c pre) [6\ 



Second time and later : 



98 

In the case of second time and later, two cases 
indicated in Figs.28A and 28B are considered. In the 
determining method of the new frequency fl^jf+ll in the 
second time and later, the frequency which satisfies the 
5 attenuation quantity dBc thus designated is preserved 
invariably within (0 The case illustrated ±n 

Fig.28A, the characteristic shown in Fig.28B is 

m {pre) \t\(D {cur) \t\ 
considered to the frequencies ^ p \}}w p L L J- 

The case illustrated in Fig.28A is one in which one 
10 side satisfies the designated attenuation quantity. In 
this case, the processing is moved to next step 41 in 
such a way as the following equation (57) . 

Further, since the frequency which satisfies the 
designated attenuation quantity always conforms to 
15 0) { / re) [t], the case that 0^ re) [f] does not satisfy does not 

exist . 

co [ ; ttr) [t + 1] = 0.5 x (cu^'H+ a^B ( 57 ) 

dB [ r ] \t + \] = dB { c pre) [t] 

The case illustrated in Fig.28B is one in which the 
20 both satisfies the designated attenuation quantity. In 
this case, the processing is moved to next step 41 in 
such a way as the following equation (58) . 
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^[r + lj-^tf] 

a>*r r) [t + 1 1 - c»jr J [r ]+ 0.5 x (^r» H- [j ) < 5 s > 

Step 45 

The coefficients of the linear-phase FIR filter is 
found from the approximated amplitude characteristic (F 
105) . 

Namely, the coefficients of the linear-phase FIR 
filter is found from the amplitude characteristic 
obtained finally. 

Fig. 29 is a view illustrating frequency response 
characteristic of the low-pass filter obtained from the 
algorithms for finding "filter that satisfies the 
attenuation quantity of the stop band and that causes the 
signal to pass the frequency point of the transition 
band" on the basis of the algorithms for finding "filter 
having the minimum starting frequency of the stop band 
while taking the filter of the zero point for avoiding 
the chessboard distortion. 

Fundamental algorithms, variables intended to find, 
and band of this case are as follows: 
Fundamental algorithms 
It is the Remez Exchange algorithms for finding the 
minimum starting frequency of the stop band that 
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satisfies the attenuation quantity of the stop band. 

• passing through (0.4 n, 12 dB) 

• 2 4 taps 

• even symmetry 

- U = 3 (filter coefficient is executed U times 
such that DC gain becomes U) 

• attenuation quantity of the stop band is not more 
than - 4 0 dB 

Variables intended to find 
It is the end frequency wp of the pass band and the 
starting frequency cos of the stop band. 

[Table 6] 

Band 



band 


Frequency region 


gain 


weight 


Pass band 


0 ^ co ^ co p 


1 


1 


Stop band 


co 3 <co^n 


0 


1 



Further, in Fig. 2 9 , the solid line satisfies the 
attenuation quantity of the stop band and indicates 
frequency response of low-pass filter causing the 
frequency point of the transition band to pass through. 

Furthermore, the black circle indicates the 
designated frequency point. Also, dotted line indicates 
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frequency which becomes zero point for avoiding the 
chessboard distortion. 

As known from Fig. 29, the low-pass filter which is 
obtained from the algorithms for finding "filter that 
causes the signal to pass the frequency point of the 
transition band and that satisfies the attenuation 
quantity of the stop band" on the basis of "filter having 
the minimum starting frequency of the stop band" has good 
frequency response characteristic. 

Subsequently, as for an algorithms on the basis of 
the algorithms for finding filter having the maximum end 
frequency of the pass band that satisfies the attenuation 
quantity of the stop band, explanation is made in 
relation to Fig. 24, Fig. 30 to Figs.33A and 33B. 

The flowchart of this case is equivalent to Fig. 24 
which is referred to in the algorithms for finding filter 
that causes the signal to pass the frequency point of the 
transition band and that satisfies the attenuation 
quantity of the stop band. 

Namely, the present algorithms utilizes the 
algorithms for finding "filter having the maximum end 
frequency of the pass band that satisfies the attenuation 
quantity of the stop band" . 

The concrete course of the algorithms, as illustrated 
in Fig. 30, cop is obtained with aforementioned algorithms 
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as inner loop, further loop for finding cos is put on 



outside thereof to determine parameter of the band. 



Namely, the outside loop causes the starting 



frequency cos of the stop band to be fixed, while the 



5 inner loop causes the end frequency cop of the pass band 
to be found that satisfies the attenuation quantity dBs 
of the stop band. 

When the amplitude characteristic thus found does not 
pass through the designated point (coc, dBc) of the 
10 transition band, cos is searched while using dichotomizing 
method like the case of searching "filter having the 
minimum starting frequency of the stop band that 
satisfies the attenuation quantity of the stop band" . 



15 however, since the algorithms for finding "filter having 
the maximum end frequency of the pass band that satisfies 
the attenuation quantity of the stop band" to be the 
inner loop is the same as that described above, thus 
description thereof is omitted. 

20 Also, the content of Fig. 24 and each step processing 

F 102, F 103, F 104, and F 105 of the algorithms 
explained below, like the case of the first approach, is 
the same as the Remez Exchange algorithms taking into 
account the frequency response of the pre-filter 

25 explained in relation to Fig . 6 . Further, the content of 



Hereinafter, concrete algorithms is indicated, 
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each processing F 206, F 207, F 208, like the case of the 
second approach, is the same as processing explained in 
relation to Fig. 11 namely the algorithms of "filter 
having the minimum starting frequency of the stop band 
that satisfies the attenuation quantity of the stop 
band". Accordingly, as for this processing, the same 
marks as Fig. 6 and Fig. 11 are used. 
Step 50 

As illustrated in Fig. 24, firstly, initial setting is 
executed (F 301) . In this initial setting, setting of the 
linear-phase FIR filter, setting of the band, setting of 
the coefficient of the pre-filter, setting of the initial 
extreme value point, designation of the attenuation 
quantity of the stop band, and setting of the initial 
frequency of dichotomizing method are executed. 

Concrete items thus set are as follows : 

• the number of tap 

• linear-phase FIR filter with even symmetry or odd 
symmetry 

• the number of the band is two 

• the starting frequency of the pass band co = 0 

• gain of the pass band 

• the end frequency of the stop band « = n 

• gain of the stop band 

• weighting to the pass band and the stop band 
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• coefficient of the pre-filter 

• attenuation quantity dBs of the stop band (namely, 
pointing largeness 52 of ripple of the stop band) 

• frequency ac of the transition band and attenuation 
5 quantity dBc thereof 

• frequency co (0) = co k <0> (k = 0, , R) that becomes 
the extreme value in the approximation band. Provided 
letter of right shoulder (i) expresses the number of 
repetition . 

10 Also, Fig. 31 is a view illustrating initial frequency 

of the dichotomizing method in the algorithms for finding 
filter that causes the signal to pass through the 
frequency point of the transition band, and that 
satisfies the attenuation quantity of the stop band. 

15 As illustrated in Fig. 31, in the present example, 

following values are given as the initial frequencies of 
the dichotomizing method. 

tt^ c "'^[o]= the vicinity of end frequency of the block 
20 band (59) 



The other side of Cuf '^[O] is one which is found with 
the algorithms for searching "filter having the minimum 
starting frequency of the stop band that satisfies the 
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attenuation quantity of the stop band" as the value 
described later. Namely, cos which is found by the above 

algorithms is taken to beC^ p '^[o]. 

CO p } ' 1^] — the vicinity of starting frequency of the 
pass band (60) 

It should be noted that cop should be made to search 
as later described equation (61) essentially, however, cop 
that satisfies the attenuation quantity of the stop band 
can not be found. Accordingly, since the minimum value of 
the cop is degree of 0.01, the starting frequency cos of 
the stop band that satisfies the attenuation quantity of 
the stop band at this time is taken to be the initial 
frequency . 

^M^qJ = the vicinity of WG (61) 

Portion of [t] at the back of cos expresses the number 
20 of cycle . 

Here, the Remez Exchange algorithms taking into 
account the frequency response of the pre-filter is 
executed to the frequency 0>j pr ^[o] , thus the maximum end 
frequency Ct^ pre *[o] of the pass band that satisfies the 
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designated attenuation quantity dBs in the stop band is 
found, the gain c/£?^ cur ^[()] corresponding to wc at this 
time is found. The following steps are explained on the 
basis of this assumption. 
5 Step 51 

The algorithms for finding "filter having the maximum 
end frequency of the pass band that satisfies the 
attenuation quantity of the stop band" is executed (F 
102, F 103, F 104, F206, F 207, F 208). 

10 Concretely, in the processing F 102, interpolation 

polynomial expression is generated for interpolating an 
amplitude characteristic from the extreme value point at 
the time of a>j C£ "^[£ ] . 

Subsequently, in the processing F 103, the new 

15 extreme value point is determined from the amplitude 

characteristic obtained from the interpolation polynomial 
expression . 

Then, in the processing F 104, repetition judgement 
of the Remez Exchange algorithms is executed. 
20 Next, in the processing F 206, the minimum 

attenuation quantity is found in the stop band. 

Next, in the processing F 207, end condition of the 
algorithms for finding filter having the minimum starting 
frequency cos of the stop band that satisfies the 
25 attenuation quantity of the stop band is obtained. 
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Also, in the processing F 208, setting of the 
starting frequency CO ^ ut ^ of the new stop band is 
executed. 

Step 52 

5 Next, the attenuation quantity of the frequency that 

is designated as the transition band (F 309) is examined. 

In the processing F 102, the attenuation quantity 
d^ tft/r )[YJ to the frequency coc which is designated as the 
transition band while using LaGrange' s interpolation 
10 polynomial expression found finally with the algorithms 
for finding "the maximum end frequency of the stop band 
that satisfies the attenuation quantity of the stop band" 
is examined. 

Step 53 

15 Comparison to the designated attenuation quantity of 

the transition band is executed (F 310) . 

Concretely, in the case that result of the comparison 
satisfies following equations (62) or (63) processing is 
moved to processing of Step 55 (F 105) . While when result 
20 of the comparison does not satisfy equations (62) or 
(63) , processing is moved to processing of Step 54 (F 
311) . 

\d^ n WdB c \{8 x (62) 

25 
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(63) 



Provided that el and e2 are extremely small values. 
Step 44 

5 Setting of the band is changed when result of the 

comparison does not satisfy equations (62) or (63) , in 
the comparison between the frequency coc of the transition 
band and the attenuation quantity dBc (F 311) . 

Concretely, the end frequency w { s ctjr) [t + l] of new pass 
10 band is made to set. Explanation is performed while 

dividing into two cases of the first loop and the second 
loop as the setting method. 
First time: 

In the case of the first time, there are three cases 
15 illustrated in Figs . 32A to 32C. 

Namely, cases illustrated in Fig.32A, Fig.32B, and 

Fig . 32C are considered to frequencies CO 

The case illustrated in Fig.32A is one in which the 
both satisfy the designated attenuation quantity dBc. In 
20 this case, the large frequency a;j C£/ ^[o] is taken to be 
solution, thus advancing toward the processing of Step 
55 . 

The case illustrated in Fig.32B is one in which the 
both do not satisfy the designated attenuation quantity. 
25 In this case, the present number of taps can not realize 
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the designated attenuation quantity. Accordingly, that 
effect is indicated to end. 

The case illustrated in Fig.32C is one in which only 
one side satisfies the designated attenuation quantity. 
5 In this case, processing is moved to Step 51 in such a 
way as following equation (64) . 

Further, in the present method for giving initial 
frequency point, the case in which <wj cl "^[o] is satisfied 
and &>^ r ^[o] is not satisfied does not exist. 

10 

wj /,rc {l]-o; a (c,,r) [o] 

^"'TO = 0-5 x (c^ c,,r) [0]+o^ re) [0]) (64) 

Second time and later: 
In the case of second time and later, two cases 
15 indicated in Figs.33A and 33B are considered. In the 

determining method of the new frequency ey^^f 4- l] in the 
second time and later, the frequency which satisfies the 
attenuation quantity dBc thus designated is preserved 
invariably within tt>j^|/4-l]. The case illustrated in 
20 Fig.33A, Fig.33B is considered to the frequencies 

co [ s pre) [t]w { s cur) [t]. 

The case illustrated in Fig.33A is one in which one 
side satisfies the designated attenuation quantity. In 
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this case, the processing is moved to next step 51 in 
such a way as the following equation (65) . 

Further, since the frequency that satisfies the 
designated attenuation quantity always conforms to 
5 &>j Ci/r) [t] / the case that W s (pr ^[f] does not satisfy does not 
exist. 

w ( J" c) [( +l] = rwj /w) [f] 

(V^" r] [t + \ | - 0.5 x {a> { r ] [* ]+ co { s pre) [f J ( 65 ) 

dB^i + l] = dB^* ) [t] 

10 The case illustrated in Fig.33B is a case that both 

satisfies designated attenuation quantity. In this case 
processing is moved to next step 51 in accordance with 
following equation (66) . 

a^'''^ (66) 
15 dB^t + l] = dB^ r) [t] 

Step 55 

The coefficients of the linear-phase FIR filter is 
found from the approximated amplitude characteristic (F 
105) . 

20 Namely, the coefficients of the linear-phase FIR 

filter is found from the amplitude characteristic 
obtained finally. 
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Fig. 34 is a diagram illustrating frequency response 
characteristic of the low-pass filter obtained from the 
algorithms for finding "filter that satisfies the 
attenuation quantity of the stop band and that enables 
5 the signal to be passed through the frequency point of 
the transition band" 

on the basis of the algorithms for finding "filter having 
:S3& the maximum end frequency of the pass band" with the 

filter of zero point for avoiding the chessboard 

: £ 10 distortion to be the pre-filter. 

Q 

m Fundamental algorithms, variables intended to find, 

and band of this case are as follows : 

: 3 

£J Fundamental algorithms 

^ It is the Remez Exchange algorithms for finding the 

^ 15 minimum starting frequency of the stop band that 

satisfies the attenuation quantity of the stop band. 

• passing through (0.4 n, 12 dB) 

• 2 4 taps 

• even symmetry 

20 • U = 3 (filter coefficient is executed U times 

such that DC gain becomes U) 

• attenuation quantity of the stop band is not more 
than - 40 dB 

Variables intended to find 
25 It is the end frequency cop of the pass band and the 
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starting frequency cos of the stop band. 

[Table 7] 

Band 

5 



band 


Frequency region 


gain 


weight 


Pass band 




1 


1 


Stop band 


a> 3 <G3^n 


0 


1 



10 Further, in Fig. 34 , the solid line satisfies the 

attenuation quantity of the stop band and indicates 
frequency response of low-pass filter causing the signal 
to be passed through the frequency point of the 
transition band. 

15 Furthermore, the black circle indicates the 

designated frequency point. Also, dotted line indicates 
frequency which becomes zero point for avoiding the 
chessboard distortion. 

As known from Fig. 34, the low-pass filter which is 

20 obtained from the algorithms for finding "filter that 
causes the signal to be passed through the frequency 
point of the transition band and that satisfies the 
attenuation quantity of the stop band" on the basis of 
"filter having the minimum starting frequency "of the stop 

25 band" has good frequency response characteristic. 
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Subsequently , the algorithms for designing of filter 
in order to find filter with the minimum number of tap 
that satisfies attenuation quantity of the stop band will 
be explained against the Remez Exchange algorithms with 
5 changeable band. 

Here, the algorithms for obtaining filter having the 
minimum number of tap that satisfies attenuation quantity 
of the designated stop band will be explained against the 
algorithms for obtaining a filter that has the maximum 
10 end frequency cop of the pass band, which end frequency 
satisfies attenuation quantity of the stop band and 
against the algorithms for obtaining a filter having the 
minimum starting frequency cos of the stop band that 
satisfies attenuation quantity of the stop band. 
15 Fig. 35 is a view illustrating flowchart of the 

algorithms for finding filter with the minimum number of 
tap that realizes attenuation quantity of the stop band. 

Here, free parameter of this algorithms, object, and 
principle of the algorithms thus enumerated are as 
20 follows : 

* free parameter: 

• it is the number of taps 

• there exists two kinds of approaches. 
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Firstly, the end frequency cop of the pass band is 
made to be variable, and the starting frequency cos of the 
stop band is made to fix. 

Secondly, the end frequency cop of the pass band 
5 is made to fix, and the starting frequency cos of the stop 
band is made to be variable . 

* object: it is an object to obtain filter having the 
minimum number of tap that satisfies attenuation quantity 
!2J dBs of the designated stop band to the algorithms in 

i 10 which one of the variables of the band is variable. 

;;n * principle: When being given no-solution in one time 

of loop, it is impossible to satisfy the attenuation 
LP quantity of the designated stop band because the number 

of tap is short, attempt is executed while increasing the 
:==a 15 number of tap by one. 

Also, the content of Fig. 35 and each step processing 
F 102, F 103, F 104, and F 105 of the algorithms 
explained below, like the case of the first approach, is 
the same as the Remez Exchange algorithms taking into 
20 account the frequency response of the pre-f ilter 

explained in relation to Fig. 6. Further, the content of 
each processing F 201, F 206, F 207, F 208, like the case 
of the second approach, is the same as processing 
explained in relation to Fig. 11 namely the algorithms of 
25 "filter having the maximum end frequency of the pass band 
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that satisfies the attenuation quantity of the stop band" 
and the algorithms of "filter having the minimum starting 
frequency of the stop band that satisfies the attenuation 
quantity of the stop band". Accordingly, as for this 
5 processing, the same marks as Fig. 6 and Fig. 11 are used. 
Step 60 

As illustrated in Fig. 35, firstly, initial setting is 
executed (F 201) . In this initial setting, setting of 
linear-phase FIR filter, setting of band, setting of 
10 coefficient of the pre-filter, setting of the initial 

extreme value point, designation of attenuation quantity 
of the stop band, and setting of initial frequency of 
dichotomizing method are executed. 

Concrete items thus set are as follows: 
15 • the number of tap 

• linear-phase FIR filter with even symmetry or odd 
symmetry 

• the number of the band is two 

• the starting frequency of the pass band oo = 0 

20 • end frequency cop of the pass band in the case of 

the algorithms for finding "filter having the minimum 
starting frequency of the stop band that satisfies the 
attenuation quantity of the stop band. 

• gain of the pass band 

25 • the end frequency of the stop band co = n 
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/ • starting frequency cos of the stop band in the case 
of the algorithms for finding "filter having the maximum 
end frequency of the pass band that satisfies the 
attenuation quantity of the stop band. 
5 • gain of the stop band 

• weighting in connection with the pass band and the 
stop band 

• coefficient of the pre-filter 

• attenuation quantity dBs of the stop band (namely, 
10 pointing largeness 52 of ripple of the stop band) 

• frequency co <0) = co k (0, (k = 0 , . . , R) that becomes 
the extreme value in the approximation band. Provided 
letter of right shoulder (i) expresses the number of 
repetition . 

15 • input of initial frequency of dichotomizing method 

Step 61 

The Reruez Exchange algorithms taking into account the 
frequency response of the pre-filter is executed (F 102, 
F 103 , F 104) . 
20 Concretely, in the processing F 102, the 

interpolation polynomial expression is generated for 
interpolating an amplitude characteristic from the 
extreme value point . 

Subsequently, in the processing F 103, the new 
25 extreme value point is determined from the amplitude 




117 

characteristic obtained from the interpolation polynomial 
expression. 

Then, in the processing F 104, repetition judgement 
of the Remez Exchange algorithms is executed. 
5 Step 62 

Next, the minimum attenuation quantity (the maximum 
weighted approximation error) of the stop band is 
obtained (F 206) . 
Step 63 

10 Discrimination is executed whether or not the end 

condition of the search of the frequency that satisfies 
the attenuation quantity of the designated stop band 
holds (F 207) . 

When the end condition holds, processing is moved to 
15 Step 67 (F 105) , while when the end condition does not 
hold, processing is moved to Step 64 (F 208) . 
Step 64 

When the end condition does not hold in connection 
with the search algorithms of the frequency that 
20 satisfies the attenuation quantity of the designated stop 
band, setting of the band is changed (F 208) . 
Step 65 

Next, comparison to the attenuation quantity of the 
designated stop band is executed (F 414) . 
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In the setting change of the band in the processing F 
208, in the case of no solution at the one time loop 
processing is moved to processing of Step 66 (F 414) . 
While with the exception of the aforementioned case, 
5 processing is returned to processing of Step 61. 
Step 66 

One tap is increased (F 415) . 

One tap is increased in addition to the present 
number of tap, thus processing is moved to initial 
10 setting processing of Step 60 (F 201) . 
Step 67 

The coefficients of the linear-phase FIR filter is 
found from the approximated amplitude characteristic (F 
105) . 

15 Fig. 36 is a view illustrating frequency response 

characteristic of the low-pass filter obtained from the 
algorithms for finding "the minimum number of tap that 
satisfies attenuation quantity of the stop band" while 
taking the filter of the zero point for avoiding the 
20 chessboard distortion. 

Fundamental algorithms, variable intended to find, 
and band of this case are as follows : 
Fundamental algorithms 
It is the Remez Exchange algorithms taking into 
25 account frequency characteristic of the pre-filter for 



10 
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finding the filter which has the minimum starting- 
frequency of the stop band, that satisfies the attenuation 
quantity of the stop band. 

• even symmetry 

• U = 3 (filter coefficient is executed U times 
such that DC gain becomes U) 

• attenuation quantity of the stop band is not more 
than - 60 dB 

Variable intended to find 

• it is the number N of tap. 

• the starting frequency cos of the stop band. 



15 



Band 



[Table 8] 
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co s <w^n 
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1 



20 Further, in Fig. 36 , the solid line indicates 

frequency response of low-pass filter of the minimum 
number of tap (14 taps) in which the attenuation quantity 
of the stop band becomes not more than —60 dB . 
Furthermore, the dotted line indicated by (D indicates 

25 frequency which becomes zero point for avoiding the 
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chessboard distortion, and the dotted line indicated by 
(D indicates period of the band given beforehand. 

As seer; from Fig. 36, the low-pass filter which is 
obtained from the algorithms for finding "the minimum 
5 number of tap which realizes attenuation quantity of the 
stop band" has good frequency response characteristic. 

Furthermore, Fig. 37 is a view illustrating frequency 
response characteristic of the low-pass filter obtained 
^ from the algorithms for finding "the minimum number of 

2 10 tap that satisfies attenuation quantity of the designated 

i=f* stop band" while taking the filter of the zero point for 

avoiding the chessboard distortion. 
iXt Fundamental algorithms, variable intended to find, 

and band of this case are as follows: 
^ 15 Fundamental algorithms 

It is the Remez Exchange algorithms taking into 
account frequency characteristic of the pre-filter for 
finding the filter which has the maximum end frequency of 
the pass band that satisfies the attenuation quantity of 
20 the stop band. 

• even symmetry 

• U = 3 (filter coefficient is executed U times 
such that DC gain becomes U) 

• attenuation quantity of the stop band is not more 
25 . than - 60 dB 
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Variable intended to find 

• it is the number N of tap. 

* the starting frequency cos of the stop band. 
5 [Table 9] 



band 


Frequency region 


gain 


weight 


Pass band 




1 


1 


Stop band 


0 . Sn^co^n 


0 


1 



Further, in Fig. 37, the solid line indicates 
frequency response of low-pass filter of the minimum 
number of tap (15 taps) in which the attenuation quantity 
of the stop band becomes not more than —60 dB . 

15 Furthermore, the dotted line indicated by indicates 
frequency which becomes zero point for avoiding the 
chessboard distortion, and the dotted line indicated by 
(2) indicates period of the band given beforehand. 

As known from Fig. 37, the low-pass filter which is 

20 obtained from the algorithms for finding "the minimum 

number of tap which realizes attenuation quantity of the 
designated stop band" has good frequency response 
characteristic . 

Subsequently, as the second modified example, 

25 explanation is performed about algorithms for finding 
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filter with the minimum numbers of tap that enables the 
signal to be pass through frequency point of transition 
band that satisfies attenuation quantity of the stop 
band . 

5 Here, algorithms for finding filter with the minimum 

number of tap that satisfies the attenuation quantity of 
the block band and and that enables the signal to be 
passed through frequency point of transition band and to 
the algorithms for finding "filter that satisfies the 
10 attenuation quantity of the block band and and that 

enables the signal to be passed through the frequency 
point of the transition band" . 

Fig. 38 is a view illustrating flowchart of the 
algorithms for finding filter with the minimum number of 
15 tap that satisfies the attenuation quantity of the stop 

band and that enables the signal to be passed through the 
frequency band of the transition band. 

Here , free parameter of this algorithms, object, and 
principle of the algorithms thus enumerated are as 
20 follows : 

* free parameter: 

• it is the number o f tap s 

• the end frequency cop of the pass band. 

• the starting frequency cos of the stop band 



123 

* object: it is an object to obtain filter having the 
minimum number of tap that satisfies attenuation quantity 
dBs of the designated stop band and that the attenuation 
quantity dBc is passed at the frequency coc of the 

5 transition band. Namely, the maximum end frequency cop of 
the pass band so as to become the attenuation quantity 
dBc at the designated frequency coc of the transition band 
and the minimum starting frequency cos of the stop band 
are determined to obtain the filter with the minimum 
10 number of tap . , 

* principle: When being given no-solution in one time 
of loop, it is impossible to satisfy the attenuation 
quantity because the number of tap is short, attempt is 
executed while increasing the number of tap by one . 

15 Further, in cases where it is impossible to realize the 
frequency point of the transition band, attempt is 
executed while increasing the number of tap by one . 

Also, the content of Fig. 38 and each step processing 
F 103, F 104, and F 105 of the algorithms explained 

20 below, like the case of the first approach, is the same 
as the Remez Exchange algorithms taking into account the 
frequency response of the pre-filter explained in 
relation to Fig . 6 . Further, the content of each 
processing F 201, F 206, F 207, F 208, like the case of 

25 the second approach, is the same as processing explained 
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in relation to Fig. 11 namely the algorithms of "filter 
having the maximum end frequency of the pass band that 
satisfies the attenuation quantity of the stop band" or 
the algorithms of "filter having the minimum starting 
5 frequency of the stop band that satisfies the attenuation 
quantity of the stop band". Further, the contents of the 
processing F 301, F 309, F 310, and F 311 are the same as 
the processing explained in relation to Fig. 24, namely, 
the same as the algorithms for finding the filter that 

10 satisfies the attenuation quantity, and, in which 

designated attenuation quantity is passed through at the 
designated frequency of the transition band. Accordingly, 
as for this processing, the same marks as Fig . 6 Fig. 11 
and Fig. 2 4 are used. 

15 Step 70 

As illustrated in Fig. 38, firstly, initial setting is 
executed (F 301) . Concretely, initial setting of the 
algorithms is executed for finding the filter that 
satisfies the attenuation quantity of the stop band, , and 

20 in which designated attenuation quantity is passed 
through of the transition band. 

Concrete items thus set are as follows : 

• the number of tap 

• linear-phase FIR filter with even symmetry or odd 
25 symmetry 
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• the number of the band ±s two 

• the starting frequency of the pass band co = 0 

• gain of the pass band 

• the end frequency of the stop band co = n 
5 • gain of the stop band 

• weighting to the pass band and the stop band 

• coefficient of the pre-filter 

• attenuation quantity dBs of the stop band (namely, 
pointing largeness 52 of ripple of the stop band) 

10 • frequency coc of the transition band and attenuation 

quantity dBc thereof 

• frequency co <0> = co k <0> (k = 0, . . , R) that becomes 
the extreme value in the approximation band. Provided 
letter of right shoulder (i) expresses the number of 

15 repetition. 

• input of initial frequency of dichotomizing method. 
Step 71 

The Remez Exchange algorithms taking into account the 
frequency response of the pre-filter is executed (F 102, 
20 F 103, F 104) . 

Concretely, in the processing F 102, interpolation 
polynomial expression is generated for interpolating an 
amplitude characteristic from the extreme value point. 
Subsequently, in the processing F 103, the new 
25 extreme value point is determined from the amplitude 
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characteristic obtained from the interpolation polynomial 
expression . 

Then, in the processing F 104 , repetition judgement 
of the Remez Exchange algorithms is executed. 
5 Step 72 

Next, the minimum attenuation quantity (the maximum 
weighted approximation error) of the stop band is 
obtained (F 20 6) . 
Step 73 

10 Discrimination is executed whether or not the end 

condition of the search of the frequency that satisfies 
the attenuation quantity of the designated stop band 
holds (F 207) . 

When the end condition holds, processing is moved to 
15 Step 77 (F 309), while when the end condition does not 
hold, processing is moved to Step 74 (F 208) . 
Step 74 

When the end condition does not hold in connection 
with the search algorithms of the frequency that 
20 satisfies the attenuation quantity of the designated stop 
band, setting of the band is changed (F 208) . 
Step 75 

Next, comparison between the attenuation quantity of 
the designated stop band and obtained attenuation 
25 quantity is executed (F 414) . 
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In change of the setting of the band in the 
processing F 208, in the case of no solution at the one 
time of the loop, processing is moved to processing of 
Step 7 6 (F 415) . In the cases with the exception of the 
5 aforementioned case, processing is returned to processing 
of Step 71. 

Step 76 

Q One tap is increased (F 415) . 

One tap is increased in addition to the present 
;S 10 number of tap, thus processing is moved to initial 

setting processing of Step 70 (F 301) . 
Step 77 

3 When the end condition of the search algorithms of 

3 the frequency that satisfies designated attenuation 

15 quantity of the stop band holds, attenuation quantity of 
the frequency which is designated to the transition band 
is examined (F 309) . 

Comparison between the attenuation quantity of the 
designated stop band and obtained attenuation quantity is 
20 performed (F 414) . 

When the end condition holds, processing is moved 
to Step 82 (F 105) , while when the end condition does not 
hold, processing is moved to Step 79 (F 311) . 
Step 79 
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When the end condition does not hold, setting of the 
band is changed (F 311) . 
Step 80 

Next, comparison between the attenuation quantity of 
5 the designated stop band and obtained attenuation 
quantity is executed (F 416) . 

In change of the setting of the band in the 
processing F 311, in the case of no solution at the one 
time of the loop, processing is moved to processing of 
10 Step 81 (F 417) . In the cases with the exception of the 

aforementioned case, processing is returned to processing 
of Step 71. 

Step 81 

One tap is increased (F 417) . 
15 One tap is increased in addition to the present 

number of tap, thus processing is moved to initial 
setting processing of Step 70 (F 301) . 
Step 83 

The coefficients of the linear-phase FIR filter is 
20 found from the approximated amplitude characteristic (F 
105) . 

Fig. 39 is a diagram illustrating frequency response 
characteristic of the low-pass filter obtained from the 
algorithms for finding "filter with the minimum number of 
25 taps which filter satisfies the attenuation quantity of 
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the stop band, and which filter enables the signal to be 
passed through the frequency point of the transition 
band" with the filter of zero point for avoiding the 
chessboard distortion to be the pre-filter. 
5 Fundamental algorithms , variable intended to find, 

and band of this case are as follows : 
Fundamental algorithms 
It is the Remez Exchange algorithms for obtaining a 
;g filter that enables the signal to be passed through the 

10 frequency point of the transition band on the basis of 
the algorithms for obtaining the filter that satisfies 
the attenuation quantity of the stop band that has the 
maximum end frequency of the pass band. 

• passing through (0.4 n, 12 dB) 
15 • even symmetry 

• U = 3 (filter coefficient is executed U times 
such that DC gain becomes U) 

• attenuation quantity of the stop band is not more 
than - 60 dB 

20 Variable intended to find 

• it is the number of taps N. 

• the end frequency cop of the pass band. 

• the starting frequency us of the stop band. 



25 
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[Table 10] 
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Further, in Fig. 39, the solid line indicates 
frequency response of low-pass filter of the minimum 
number of tap (17 taps) in which the attenuation quantity 

10 of the stop band becomes not more than — 60dB , and the 

attenuation quantity becomes not more than — 12dB on the 
frequency 0 . 4n of the transition band. Further, dotted 
line indicates the frequency that becomes zero point in 
order to avoid the chessboard distortion, the black 

15 circle indicates the designated frequency point of the 
transition band . 

As seen from Fig. 39, the low-pass filter which is 
obtained from the algorithms for finding "filter with the 
minimum number of taps which filter satisfies attenuation 

20 quantity of the stop band, and which filter enables the 
signal to be passed through the frequency point of the 
transition band" has good frequency response 
characteristic . 

As described-above, according to the present 

25 invention, there is an advantage that is possible to 
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stably preserve weighted approximation error in such a 
way that ecpai-ripple is not broken in comparison with the 
conventional method. Further, there is an advantage that 
is possible to preserve gain of the pass band into a 
constant value in comparison with the conventional 
method. 

While the invention has been described with reference 
to specific embodiment chosen for purpose of 
illustration, it should be apparent that numerous 
modifications could be made thereto by those skilled in 
the art without departing from the basic concept and 
scope of the invention. 



